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The di-electron width of an Upsilon meson is the decay rate of the Upsilon into an 
electron-positron pair, expressed in units of energy. We measure the di-electron 
width of the Upsilon(lS) meson to be 1.354 ± 0.004 ± 0.020 keV (the first un- 
certainty is statistical and the second is systematic), the di-electron width of the 
Upsilon(2S) to be 0.619 ± 0.004 ± 0.010 keV and that of the Upsilon(3S) to be 
0.446 ± 0.004 ± 0.007 keV. We determine these values with better than 2% pre- 
cision by integrating the Upsilon production cross-section from electron-positron 
collisions over their collision energy. Our incident electrons and positrons were 
accelerated and collided in the Cornell Electron Storage Ring, and the Upsilon 
decay products were observed by the CLEO-IH detector. The di-electron widths 
probe the wavefunctions of the Strongly-interacting bottom quarks that consti- 
tute the three Upsilon mesons, information which is especially interesting to check 
high-precision Lattice QCD calculations of the nuclear Strong force. 
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Chapter 1 

Introduction and Motivation 

1.1 The Y Di-electron Width and Why it is Important 

An Upsilon (T) meson is a composite particle consisting of a bottom quark (6) 
and an anti-bottom quark (b) bound with their spins ahgned in a J = 1 quantum 
mechanical wavefunction, where J is the total angular momentum. This meson is 
a nuclear analogy of ortho-positronium in atomic physics. The di-electron width 
is the rate of T decay into an electron/positron pair, and measuring it provides 
unique experimental access to the physical size of the bb wavefunction and its total 
decay rate — the average extension of the T meson in both space and time. 

The bb system, also known as bottomonium, is the most non-relativistic system 
of quarks bound by the nuclear Strong force. This is because the bottom quark 
is the heaviest quark that can participate in the Strong force, the top decaying 
immediately into bottom by the Electroweak force. Unlike much more abundant 
protons and neutrons, whose masses consist almost entirely of the kinetic energy 
of the constituent quarks and gluons, 94% of the mass of the lightest T consists of 
the mass of its two bottom quarks. This simplifies the dynamics of bottomonium 
and even permits description in terms of a potential, making it a good testing 
ground for Strong force calculations. 

Quantum Chromodynamics (QCD) has long been accepted as the correct de- 
scription of the nuclear Strong force (with possible modifications only at TeV en- 
ergies and above) because of its success in predicting scattering interactions above 
one GeV and its qualitative explanation of low-energy phenomena like quark con- 
finement. Today, the Lattice QCD technique, which simulates QCD on a computer. 



is yielding few-percent calculations of low-energy phenomena from first principles — 
in particular, T properties such as the di-electron widths. Precise experimental 
knowledge of the T di-electron widths will test the new Lattice QCD techniques 
that made this advance possible. 

The di-electron widths check Lattice QCD in a way that is key for Electroweak 
physics. The CP violation parameters Vtd and V^b, fundamental constants in the 
Standard Model, could be extracted from existing hadronic measurements much 
more precisely if the strength of the force between quarks were better known. 
Lattice QCD can help, but precise Lattice results will only be trusted if similar 
calculations can be experimentally verified. The T di-electron width closely resem- 
bles the factor that limits our knowledge of Vtd, and thus will provide a cross-check 
that will either lend credence to or cast doubt on the Vtd extraction. 

Di-electron widths of the T resonances have been measured before, but not with 
the precision that is now being demanded by Lattice QCD. This document repre- 
sents a comprehensive study of the T(IS'), T{2S), and T(3S') di-electron widths, 
with 50 times the data of any previous measurement. We present di-electron width 
measurements of the T(IS'), T(2S'), and T(3S') with 1.5%, L8%, and 1.8% total 
uncertainty, respectively. This is the second-ever measurement of the T{3S) di- 
electron width, improving its precision by a factor of five. Furthermore, measuring 
all three resonances in the same study permits us to derive very precise ratios of 
di-electron widths, where the tightest constraint on theory is likely to be. Without 
this measurement, comparisons with Lattice QCD would probably be limited by 
experiment. 



1.2 The BottomoniurQ Potential and Mass Eigenstates 

The 6-quark and the 6-quark in bottonionium attract each other by the nuclear 
Strong force, which in QCD is mediated by gluons, the nuclear analogy of virtual 
photons. The two quarks are charged with opposite "colors," in a quantum me- 
chanical superposition of red/anti-red, blue/anti-blue, and green/anti-green states, 
which are constantly traded for each other by the doubly-colored gluons. (The 
exchange of a red/anti-green gluon will turn a red/anti-red bb system green/anti- 
green, for instance. See Figure 1.1(a).) Because the gluons carry color charge, they 
can interact with other gluons and spawn complicated networks of interactions be- 
tween the two quarks (Figure 1.1(b)), which increases the interaction strength with 
distance. Bottom quarks are usually separated by a femtometer, and at this dis- 
tance scale, the coupling constant of QCD is of order unity. Feynman diagrams 
with many vertices are not suppressed relative to simple diagrams, and therefore a 
calculation of the force between the two quarks does not submit to a perturbative 
expansion. 

A very successful model of the force between b and b consists of a Coulomb-like 
potential at short distances (though about 50 times stronger) and a linear potential 
at large distances (which limits to a constant force of about 14 tons), as illustrated 
in Figure 1.2(a) [1]. At large distances, a string of self-interacting gluons, stretched 
between the two quarks, is responsible for the linear component. This string will 
generate a real quark/anti-quark pair and snap if stretched with sufficient energy. 
There are three J = 1 solutions to the Schrodinger equation below this threshold: 
they are labeled T(IS'), T(2S'), and T(3S'). States above this threshold have a 
very different pattern of decay modes and are beyond the scope of this study. 

These three mass eigenstates are the bottomonium equivalent of atomic energy 
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Figure 1.1: (a) An example of a gluon as a force propagator and a carrier of color 
charge. Directed lines are quarks and springs are gluons. (b) An exchange between 
two quarks involving a complicated network of gluons and a light quark/ ant i- quark 
pair. Grey shading indicates the sum of all such amplitudes. 
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Figure 1.2: (a) Schematic of the potential energy between two quarks as a function 
of separation, (b) Quantitative level diagram oi J = 1 bb solutions (T). The T{4S) 
lies just above the string-breaking threshold (dashed line). 



levels — discrete lines of allowed mass-energies. But, just as in atomic spectra, the 
short lifetimes of these states imply a broadening of their spectral lines: they are 
not perfect time-independent solutions. The full-width of an T resonance at half- 
maximum, r, is equal to its decay rate, in analogy with excited atomic resonances. 
If we partition the T decays into distinguishable modes, one being T — > e"'"e~, the 
total decay rate is a sum of those modes. Hence, Fee = ^ Bee, where Tee is the 
di-electron width and Bee is the fraction of T mesons that decay to e^e~, that is, 
the branching fraction to e+e^. 

The T meson decays into e+e^ by bb annihilation (Figure 1.3), which is a 
point-like interaction. The b and the b must fluctuate to the same point in space 
for the reaction to proceed. This probability, which is the square of the bb spatial 
wavefunction evaluated at the origin (|-?/'(0, 0, 0)p), is therefore a factor in Fee. 

r„ = 30.^i5^MM^ ,1.1) 

^ 3 Mr 

where Qb = 1/3, the 6-quark electric charge, a is the Electromagnetic fine structure 
constant and Mx is the T mass [2]. This is a non-relativistic approximation: 
relativistic corrections replace the wavefunction at the origin with an integral of 
values very close to the origin. Because of this dependence on knowledge of the 
bb wavefunction, and therefore the potential, a first-principles calculation Fee will 
require non-perturbative QCD. 

1.3 Lattice QCD 

Feynman path integrals provide a general approach to quantum field theory that 
don't rely on a perturbative expansion. In this formalism, the amplitude of a 
process is calculated as a weighted sum of all possible ways it can proceed. The 
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value of every field at every point in space in the initial state is allowed to vary 
as an arbitrary function of time to the final state, and these paths are weighted 
by their action. This is a generalization of Lagrange's method in classical physics, 
in which the true path is the one which minimizes action. In quantum physics, 
all paths which nearly minimize action contribute to the amplitude of a process 
(Figure 1.4(a)). 

To calculate a sum over a set of arbitrary paths, one must discretize space- 
time into time slices and space cubes. The path of a field value in a space cube 
from the initial state to the final state is a sequence of values for each time slice 
(Figure 1.4(b)). A sequence of A^ values is a vector in an A^-dimensional vector 
space: the weighting factor is integrated over these vector spaces. To obtain a 
realistic result, one must afterward limit the discretization scale to zero. 

In general, realistic problems like QCD, the integral will not be analytic and 
must be solved by numerical integration. The integral will have a fixed number of 
dimensions, which implies a fixed discretization of space-time that can only be lifted 
by extrapolating several calculations toward zero lattice size. This discretization is 
the lattice of Lattice QCD. Quark field values are represented on a four- dimensional 
lattice of space-time points with gluon field values on the edges connecting them. 

This is a very computationally intensive problem, since the number of dimen- 
sions in the integral scales with the number of grid points, and one must maximize 
the number of grid points to extrapolate to the continuum limit. Over the past 30 
years, theorists have improved the algorithms of Lattice QCD and sought approx- 
imations to make realistic calculations tractable. 

The most time-consuming part of most Lattice QCD calculations is the po- 
larization of the vacuum by light quarks. In terms of Feynman diagrams, these 
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Figure 1.3: Decay diagram of T{nS) — > e+e . The Z^ contribution is 0.25% of the 
total rate. 
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Figure 1.4: (a) Paths in field strength {ip) versus time. The minimum (white) is 
the classical solution and the path which contributes the most to the quantum am- 
plitude; shades of grey represent other quantum paths with smaller contributions 
to the total amplitude, (b) A path approximated by a discrete time sequence. 




u, d, s 



Figure 1.5: Vacuum polarization by light quarks. 

are interruptions of a gluon propagator by loops of uu, dd and ss pairs, which 
can be ignored or suppressed by assuming infinite or very heavy up, down, and 
strange quark masses (Figure 1.5). This approximation is known as the quenched 
approximation, and it permits calculations with 10-20% systematic uncertainties. 

This situation was dramatically improved in the late 1990's by the development 
of new algorithms based on the Symanzik-improved staggered-quark formalism [3]. 
These algorithms are by far the most efficient known, and the formalism features an 
exact chiral symmetry that particularly benefits simulations with small light quark 
masses. Realistic up and down quark masses are still out of reach, but simulations 
using masses three times too large can be accurately extrapolated with chiral 
perturbation theory. Thus, "unquenched" calculations are now possible, resulting 
in the accurate prediction of many masses and decay rates, as demonstrated in 
Figure 1.6. 

This algorithmic speed comes at a conceptual price: the staggered-quark for- 
malism introduces four equivalent species of each quark field, called "tastes." These 
are artifacts of the formalism and are unrelated to quark flavor. Each of these tastes 
contributes to the vacuum polarization, resulting in loop contributions which are 
four times too large. To correct for this, the quark determinant in the action is 
replaced by its fourth root, a procedure which is rigorous in the free-field theory 
and in perturbative QCD, but introduces violations of Lorentz symmetry at short 
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Figure 1.6: Lattice QCD calculations divided by experimental measurements for 
nine quantities, without and with quark vacuum polarization (left and right panels, 
respectively). (The notation T(n — 15*) refers to the mass difference between T(n) 
and T{1S).) 
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distances in the lattice simulation. Although these non-physical effects can be 
removed by interpolating between the lattice points with perturbative QCD, this 
issue makes the new algorithms controversial, and it is one aspect that will be 
tested by confrontations with experiment. 

The di-electron width may be determined from T simulations by extracting 
the bb wavefunction at the origin and applying Equation 1.1. In a path integral 
context, the wavefunction is the quark field amplitude. These simulations employ 
a non-relativistic QCD action with relativistic corrections (NRQCD), because the 
de Broglie wavelength of a massive b quark would require impractically narrow 
time slices. 

Simulations of the T mesons have been generated by the HPQCD collaboration, 
but the determination of Fee from them is not yet complete [4]. To properly 
calculate Fee, one needs to correct the lattice wavefunction for discretization effects 
through a constant, Z^^l°^, that matches the lattice approximations of the virtual 
photon current to a continuum renormalization scheme. The leading term in ^^"atch 
is on the order of the Strong coupling constant as, about 20%. This calculation 
is in progress. However, ^™atcii largely cancels in ratios of Tee'- for instance. Fee 
from the T(2S') divided by Fee from the T(IS') can already be extracted with a 
10% uncertainty. 

= 0.43 ±0.04. 1.2 



Fee(25) 



ree(15) 



HPQCD 



This uncertainty is primarily due to residual discretization errors, evident from 
the steep dependence of the result on lattice spacing size (Figure 1.7). When the 
discretization correction has been calculated, the uncertainty in this ratio should 
be only a few percent, while absolute Fee values should have uncertainties on the 
10%-level. This is why experimentally precise ratios of Fee are also valuable. The 
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ratios test the NRQCD treatment of the b quarks with high precision, though they 
are less sensitive to the ^match corrections. 



1.4 Relationship to Electroweak Parameter Extraction 

Lattice verification of Fee is particularly significant for an application of the tech- 
nique to Electroweak physics. Vertices in Feynman diagrams that join a top quark, 
a down quark, and a W boson contribute an a priori unknown factor, Vtd, to the 
amplitude (Figure 1.8). This parameter is a fundamental constant in the Stan- 
dard Model and is essential for violation of charge-parity (CP) symmetry: if Vtd 
were zero, the Standard Model would be CP symmetric (exchanging particles for 
antiparticles and mirror-transforming space would preserve all observables) . To 
determine Vtd, one must resort to measurements of bound quark systems, because 
bare quarks do not exist in nature. The transition rates for these systems depend 
both on Vtd and on QCD factors related to the structure of the bound state. Lattice 
QCD can calculate these factors and thereby extract Vtd- 

The most sensitive probe of Vtd is B^-B'^ mixing. A B^ meson is a bound 
state of d and h quarks, and a B^ meson is its charge conjugate, hd. These two 
mesons can mix, that is, a B^ can transform into a B^ and vice-versa, through the 
diagram illustrated in Figure 1.9(a). The heavy top quark dominates in this loop 
and provides a factor of Vtd for each vertex with a down quark. The rate of this 
process is extremely well-known: 0.509 ± 0.004 ps^^, a 1% measurement [5]. 

Despite this precision, the Vtd extraction has 20% uncertainties from Strong 
interactions. To illustrate the infiuence of the Strong interaction on B^-B^ mixing, 
we re-draw the Feynman diagram as a space-time diagram in Figure 1.9(b). The 
W-t loop is a very short-range process (~0.001 fm). By comparison, the average 
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Figure 1.7: HPQCD calculations of ree(2S')/ree(15') times the ratio of masses 
squared as a function of lattice grid size a. Square data points represent calcula- 
tions with light quark masses close to their natural values, and triangular points 
are quenched (infinite light quark masses). 




Figure 1.8: The vertex joining top, down, and W . 

distance between the b and d quarks is set by the QCD potential (~fm), just as it 
is for bb in the T meson. Just as in Fge, the rate of B^-B^ mixing is determined 
by the probabihty that the two quarks will fluctuate to the same point in space, 
and this probability is characterized by the B meson decay constant Jb- 



T[B ^ B) = (known factors) x 



2 



;i.3) 



fB^BB^Vt/ 
QCD 

The B^-B^ mixing amplitude depends on two factors of J'b, one from the B^ wave- 
function and the other from B^ (see Figure 1.9(b)). Another factor, known as the 
Bag parameter Bb, corrects for gluons connecting the B^ and 5" (Figure 1.9(c)). 
Its uncertainty is more easily controlled. Our knowledge of Vtd is therefore domi- 
nated by the uncertainty in j'b- 

In principle, one can measure Jb experimentally through B^ -^ fi^u or r+z/, 
illustrated in Figure 1.10. The charged B~^ has different quark content from the 
neutral -B°, but its decay rate depends on /^ because up and down quark masses 
are both much smaller than the bottom quark mass, and flavors do not enter the 
QCD calculation. Unfortunately, this process is suppressed by Vub, to the extent 
that it has yet to be observed in 88.9 million B"^ decays at BaBar [6]. Given the 
low rate of this decay and the challenge of reconstruction, it will take a long time 
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Figure 1.9: (a) One of the two diagrams dominating B^-B^ mixing (the other 
exchanges t ^^ W). (b) The same diagram emphasizing the Strong force between 
the quarks and the relative size of the t-W loop, (c) Diagrams that contribute to 
the Bag parameter, which is not a part of j'b- 
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to accumulate enough data to make a statistically significant measurement of fs- 
The B decay constant may instead be extracted from Lattice QCD simulations 
of B mesons in much the same way as Fee is from T simulations: by sampling the 
wavefunction at the origin. The HPQCD collaboration has found J'b to be 216 ± 
22 MeV (see Figure 1.11) [7]. Like the T, the B meson is modeled with NRQCD, 
and the discretization issues and corrections to this calculation are analogous to 
Fee- The largest uncertainty in fs is in the Z^^^^^^ constant that matches lattice 
approximations of the virtual W^ current to a continuum renormalization scheme. 
This ■Z^m1itei7'^'^*°'^ ^^^ been calculated to first order in ag, but uncertainty in the a^^ 
term imposes a 9% uncertainty on j'b, which dominates the 10% uncertainty cited 
above. Finer calculations of j'b are in progress. 

Lattice calculations of fs would be viewed with suspicion if calculations of Fee 
do not match experiment at a comparable level of precision. From the lattice's 
perspective, the only difference between these two calculations is the mass of one 
of the two quarks: a bottom quark is replaced by a light quark. This is not a trivial 
distinction: it is also worthwhile checking the lattice calculation of the D meson 
decay constant, in which a charm quark and a light quark annihilate, with experi- 
mental results from CLEO-c that are now becoming available [8]. The D meson is 
a heavy/light quark combination, just like the B meson, so fn is physically more 
analogous to j'b than Fee is- However, the D meson is a more relativistic system, 
the charm quark being four times lighter than bottom, so instead of simulating 
charm quarks with NRQCD, the D meson simulations use a relativistic approxi- 
mation known as the FermiLab action [3]. Thus, Fee tests the treatment of heavy 
quarks in the Jb calculation and j'd tests the heavy-light simulation and matching 
the virtual W^ current to the continuum with ^m'atch- Ratios of Fee are particu- 
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Figure 1.10: Decay diagram of B^ —>■ t^v. 
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Figure 1.11: HPQCD calculations of fsV^B as a function of the light quark mass 
used in the simulations. The solid curve extrapolates from the simulations to the 
natural rriq/ms of 1/27.4. 
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larly applicable to this test, since they will be more sensitive to the treatment of 
heavy quarks than to ^mTtch i^ ^ee- Experimental verification of Fee and Fee ratio 
calculations are therefore key to our confidence in fs and the extraction of Vtd- 



Chapter 2 
Measurement Technique 

2.1 Scans of T Resonance Production 

To determine the decay rate of T — ;► e"'"e~, we use a special strategy available to 
e^e~ colliders: we measure the total cross-section of e^e" -^ T, the time-reversed 
process. This cross-section is related to the hb wavefunction at the origin for the 
same reason as T^e (Figure 2.1, an application of crossing symmetry). 

/ a(e+e- -. T) dEcM = 3Q,' 647rV '^^°; °',°^'' (2.1) 

J Mx 

where the dEcu integration is performed over e~^e~ center-of-mass energies [2]. 

To obtain Fee in terms of /cr(e+e^ — ^ T) dEcM, we combine the above with 

Equation 1.1. 

r- = TTT / ^(e+e- ^ T) dEcu- (2.2) 

This is more general than Equations 1.1 and 2.1, which only apply in the non- 
relativistic limit. In the fully relativistic case, |-?/'(0, 0, 0)p must be replaced with 
an integral of wavefunction values near the origin, which cancels in Equation 2.2. 
This may seem like a very indirect way of measuring Fge. Why do we not count 
T — > e~^e^ decays relative to the number of T mesons produced, for instance? 
The reason is because such a fraction would be Bee, rather than the decay rate. 
We would need to multiply by F, the rest mass distribution of the T meson, to 
determine Fee, and this is experimentally inaccessible: F is on the order of 50 keV, 
which is about a hundred times narrower than the beam energy spread of an e+e^ 
collider and a thousand times narrower than detector resolution. Equation 2.2 
provides direct access to Fee, which, as a collateral benefit, can be combined with 

18 
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Bee to obtain T. 

To evaluate J o"(e+e~ -^ T) dEcu, we fit a curve to the T lineshape, that is, 
the T production cross-section as a function of e+e~ colhsion energy. We then 
integrate this curve analytically. To construct our fit function, we begin with the 
natural lineshape of the T(IS'), T{2S), and T{3S) resonances, a Breit-Wigner: 

.(e-e-^T)(BcM)^(|;r„)i^_Ig^^^^ (2,3) 

The observed spectrum is smeared by a unit Gaussian spread in incident beam en- 
ergies (~4 MeV), as discussed above. We represent this smearing by a convolution, 
but the integral is unchanged. 

The high-energy side of the lineshape is also distorted by initial-state radiation 
(ISR): e+e^ -^ T events are hard to distinguish from e+e" -^ 7T for sufficiently 
small photon energies E^. These events add to the apparent cross- sect ion, and 
contribute a high-energy tail that falls of as l/E^, causing the integral to diverge. 
We could introduce an artificial cut-off, but then the Fee we report would be a 
function of that cut-off. Instead, we include the ISR distortion in our fit function 
to match the data, but report the integral with no ISR contribution, a procedure 
depicted in Figure 2.2. This means that a{e^e^ -^ T) in Equation 2.2 represents 
the T cross-section with no ISR photons at all, and the Fee we derive is devoid of 
final-state radiation (T -^ 7e+e^). 

In addition to e^e^ -^ T and 7T, electron-positron collisions in the 9.4- 
10.4 GeV range can also undergo the following continuum processes, which also 
contribute to the observed cross-section: 

a. qq, n~^n~ , and t^t^ purely through annihilation (s-channel. Figure 2.3(a)), 

b. Bhabha e'^e~ through annihilation (s-channel) and Coulomb scattering 
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Figure 2.1: Diagrams of T ^ e^e and e+e — > T. Both feature the same b-b-'y/Z^ 
vertex whose rate is set by |-?/'(0, 0, 0)p. 
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Figure 2.2: The anatomy of an T lineshape scan (cross-section versus center-of- 
mass energy), including the natural lineshape (dashed peak), beam energy spread 
and ISR tail (solid), and backgrounds (vertical offset). 
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Figure 2.3: Survey of continuum backgrounds: (a) s-channel fermion pair pro- 
duction, // may be qq, fJ,^fi', or t^t^ , (b) t-channel exchange which dominates 
Bhabha (e^e~) scattering, (c) e~^e~ annihilation into two real photons (the other, 
undrawn diagram simply exchanges the identity of the two outgoing photons), 
(d) fusion of two virtual photons into a low-momentum hadronic state accompa- 
nied by high-energy e+e^. 
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(t-channel, Figure 2.3(b)), 

c. 77 through annihilation (Figure 2.3(c)), and 

d. e~^e^X via the fusion of two virtual photons from a grazing collision (Fig- 
ure 2.3(d)). 

At these energies, muon- and tau-pair cross-sections are 1 nb, and qq are 4 nb, 
decreasing with center-of-mass energy as 1/s (s = -Ecm^)- Bhabha events are 
by far the most abundant; in fact, the Bhabha cross-section diverges if glancing- 
angle scatters are included. The 77 cross-section diverges also, but less rapidly 
as a function of angle. Bhabha and 77 cross-sections both decrease as 1/s. The 
last process, two-photon fusion, generates low-momentum hadronic particles X 
and two electrons (e+e^), at least one of which grazes the incident beam-line. 
The two-photon fusion cross-section increases with center-of-mass energy, but very 
slowly, as logs. Some of these non-T processes can be hard to distinguish or 
are indistinguishable from T events, and therefore can be confused with signal. 
Fortunately, the continuum cross-section is a much smoother function of Ecu than 
the T cross-section, so the T peak appears to stand on a fiat continuum plateau, 
also depicted in Figure 2.2. We add terms to the fit function to accommodate 
these effects as well. 

When a continuum final state is truly indistinguishable from an T decay, as 
is the case for e~^e~ -^ qq and e~^e~ ^ T ^ qq, the cross-sections don't simply 
add. The complex amplitudes for these processes add, the square of which is 
proportional to cross-section: 

I 1 2 I 1 2 I 1 2 

'^rcs+contl,-C!'CMJ OC UAfcs T ^cont U^rcs T U^cont T ^/v6(^^i-qs ^contjj v-^'^j 
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so we can re-write o"res+cont as 

0"rcs+cont(-EcM) = O'rcs(-EcM) + 0"cont(-EcM) + 5"int(-EcM) (2.5) 

where a denotes cross-section and A amplitude, "res" for the resonant (T) con- 
tribution and "cont" for the continuum. The interference term, crint, is a function 
of -EcM, hke the famihar a^s and (Tcont, but it can be negative and sometimes 
larger than (Tj-cs- That is, introducing another way for e^e~ to produce qq can 
actually decrease the qq cross-section! We calculate this interference term from a 
Breit-Wigner amplitude (Feynman propagator) of the form 

Acs oc — -— (2.6) 

EcM-Mr + iT/2 ^ ' 

and a constant continuum with constant phase 0o (resonance phase minus contin- 
uum phase at i?cM ^ Mx). We find 

^intiEcu) = ttintO-rcs(^CM) f 2 ^^ ^ COS 00 + slu 0o J (2.7) 



where aint = y'^cont (MxVSree) Tj. (2.8) 

The magnitude of each interference correction is characterized by ai^t, which is 
a constant derived from the continuum cross-section, the resonance magnitude 
Fee, and the decay rate Fj to the given final state (in this case qq). Given the 
continuum qq cross-section (through R = a{e^e^ -^ qq)/a{e^e^ -^ /^^A^ ) [9]) 
and the resonance qq branching ratio F^^/Fee (assuming Bee = B^^, this is R), we 
calculate ajnt to be 0.0186, 0.0179, and 0.0182 for the T{1S), T{2S), and T{3S), 
respectively, with 3% uncertainties. Note that if 0o = 0, o"int < below Mr and 
aint > above Mr- If, however, 0o = =t7r/2, cTint will have exactly the same energy 
dependence as (Xres, and thus be indistinguishable from the T cross-section itself. 
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Continuum Bhabhas, fi^fi^, t^t^ , and qq all interfere with e+e^ ^ T ^ // 
with phase angle 0o = 0. We can see this by considering that all of these processes 
are purely QED except for the formation, propagation, and disintegration of the 
T resonance. The tree-level QED amplitudes are real because they all feature an 
even number of photon vertices (each of which contributes a factor of i). While 
T production introduces a factor of the conjugated wavefunction at the origin, 
'?/'(0, 0, 0)*, the T decay part of the diagram multiplies it by ■?/'(0,0,0). The T 
propagator (Equation 2.6) is real for Eqm <^ My. Therefore, 0o = or 0o = tt. 
We see in a scan of T(IS') -^ fi^fi^ (Figure 2.4) that interference is destructive 
below resonance and constructive above, which indicates 0o = 0. 

2.2 Y Final States and Hadronic Cross-section 

In our e~^e~ collisions, T mesons are produced nearly at rest and immediately 
decay. We only ever observe the T decay products. An T may decay into 

a. leptonic final states: e~^e~, fi'^fi~, and t^t^, through an s-channel virtual 
photon (or Z^, with 1.5% contribution to the rate), 

b. hadronic final states via the hadronization of qq, ggg, or gg'~f, 

c. lower-mass bb states, accompanied by pions or photons, 

d. neutrino pairs exclusively through Z^, and 

e. possibly other, exotic, modes. 

The fi^fi" branching fractions, i3^^, have been measured with 2-5% precision 
for the T(IS'), T(2S'), and T(3S') [12], and the e^e^ and r^r^ decays are expected 
to have the same amplitudes as fi^fi'^. Thus, the branching fractions. Bee, -B^/i, 
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Figure 2.4: Efficiency-corrected, background-subtracted fi^fi^ cross-section versus 
Ecu, fitted with (a) resonance, continuum, and interference normalization allowed 
to float independently, and (b) interference fixed with the opposite phase. The fit 
prefers 0o = over (pQ = n hj 18.5 standard deviations. 
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and Brr, are nearly equal, with only a tiny correction from lepton mass, which is 
0.05% for the heavy r lepton. This assumption is called Lepton Universality. 

States containing bare quarks or gluons (partons), like qq, ggg, and gg'y, must 
hadronize before propagating to the detector. That is, strings of self-interacting 
gluons, stretched between the partons, will generate new quark/anti-quark pairs 
when stretched sufficiently far. These new quarks will clothe the original partons, 
such that a macroscopic detector will only ever observe mesons {qq bound states) 
and baryons {qqq or qqq). Hadronization is a random process, leading to a broad 
spectrum of event topologies, with as many as twenty particles in the final state. 
Most T mesons decay hadronically. 

Only the T(2S') and the T(3S') have appreciable decay rates to other bb states. 
(The T(IS') -^ rib{lS) branching fraction is expected to be less than 10"*^ [13].) 
These decays are the bottomonium equivalent of atomic transitions, but in addi- 
tion to emitting monoenergetic photons in A J = 1 decays, bottomonium can emit 
monoenergetic vrvr (charged or neutral) and 77 when decaying with AJ = 0. Fig- 
ure 2.5 plots the energy levels of the most well-known bb states and the transitions 
between them. 

The Z^ boson at T masses is 80 GeV off-shell, while the photon is only 10 GeV 
off shell, so the Z° contributes to 1.5% of the T meson's Electroweak decays (e^e^, 
/x"'"/i~, T~^T~ , and qq). The Electroweak decays account for (3 -|- R)B^^ = 16% of 
all T decays, where R, the branching ratio of qq to fi'^fi^, has a value of 3.58 ± 
0.14 [9]. The T -^ Z^ branching fraction is therefore 0.25%, and T ^ Z^ ^ ui? is 
0.05%, which is negligible at our level of precision. 

Finally, we do not exclude the possibility that unknown T decays exist, or that 
their branching fractions are larger than a few percent. These modes may resemble 
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Figure 2.5: States and transitions in the bb system. Only J = 1 T mesons are produced directly by e'^e , and of those, only 
T(2S') and T{3S) decay significantly into lower bb states. 
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hadronic decays, or have exotic signatures that have been overlooked. 

For the sake of this analysis, we will classify T decay modes as "leptonic," 
meaning e"'"e~, fi^fi~, and r^r^ exclusively (no uu), or "hadronic," meaning ev- 
erything else. By this designation, there are "hadronic" final states that contain 
no hadrons at all, such as the T(2S') -^ jXbii^P) -^ 77^(15*) -^ 77e+e~ chain, 
T -^ uu, and T -^ wimp wimp (where wiMPs are cosmologically-motivated invis- 
ible particles). This classification is convenient and has been used by previous Fge 
analyses. 

Experimentally, the T cross-section is the number of e^e^ -^ T events that 
occurred divided by the time-integrated luminosity of the e~^e~ collisions. To iden- 
tify e"'"e~ -^ T events, we select events that look like hadronic T decays, because 
the leptonic final states are hard to distinguish from continuum backgrounds and 
account for only 6-7.5% of the T decays. Thus, we count e^e^ ^ T ^ hadronic 
events and our cross-section is a{e^e^ ^ T ^ hadronic). This hadronic cross- 
section is a constant multiple of the total cross-section 

a{e^e^ ^ T ^ hadronic) = a{e^e^ ^ T) x Fhad/rtot (2.9) 

for all Ecu- We fit our lineshape function to hadronic cross-section versus -Ecm 
data and thereby derive FeeFhad/rtot- To obtain T^e, we divide by i3had = Thad/rtot, 
which is (1 — Bee — B^^ — Brr) by definition. Applying Lepton Universality, we use 

T^ f eef had/f tot 



ee 



1 - 3i3^^ 



(2.10) 



to take advantage of the well-measured B^^. With Tee, we again assume Bee = -B^/i 
to calculate F = Tee/B^^. 

The upcoming chapters will each present one aspect of the FegFhad/rtot niea- 
surement. 
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Chapter 3 will describe the e^e^ collider and particle detector that were used to 
generate and count e+e^ ^ T ^ hadronic events. 

Chapter 4 will define the e^e~ ^ T ^ hadronic selection criteria and explain 
how background events are removed from that count. 

Chapter 5 will derive the correction for hadronic T events missing from the sam- 
ple, that is, the efficiency of the selection criteria. 

Chapter 6 will explain how we measure integrated luminosity, thereby converting 
our hadronic event counts into hadronic cross-sections. 

Chapter 7 will show how we use cross-section data to put an upper limit on the 
uncertainty in beam energy measurements. 

Chapter 8 will describe the fit function and fit results in detail, and 

Chapter 9 will present all Teerhad/rtot, Tge, T, and |-?/'(0, 0, 0)p results. 



Chapter 3 



Collider and Detector 

In this Chapter, we present the apparatus we used to colhde electrons and positrons 
and collect T decay products. 

3.1 Cornell Electron Storage Ring 

Our electron-positron pairs collided in the Cornell Electron Storage Ring (CESR), 
a 768 m-circumference, symmetric storage ring and collider in Ithaca, NY. This 
collider covers a very broad range of e+e~ energies, from the charmonium region at 
Eheam = -Ecm/2 = 1.8 GcV through the excited T states at 5.5 GeV. The T(IS'), 
T(2S'), and T{3S) masses, between 4.7 and 5.2 GeV, lie in CESR's optimal range. 
In fact, scans of the T(IS') through T{3S) are among the first data taken by CESR 
in 1979 (Figure 3.1). 

Copper dipole magnets confine the beams to their orbits, alternating with 
quadrapole and sextapole magnets for focusing. Superconducting quadruples pro- 
vide the final focusing of the beams, only 30 cm from the interaction point, allowing 
the collisions to reach instantaneous luminosities of 10^^ cm^^ s^^. Like all syn- 
chrotrons, the beam is pulsed to permit acceleration: beam bunches are timed to 
enter radio frequency (RF) standing waves just when the electric field is maximal. 
In CESR, nine trains of five bunches circulate in the ring at once, with 1.15x10^*^ 
particles per bunch. Both the electron beam and the positron beam are enclosed in 
the same beam-pipe, so they need to be separated electrostatically. The resulting 
orbit is called a "pretzel orbit" because the beams twine around each other like 
twisted pair cable. 
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We determine the beam energy by measuring the magnetic field in two dipole 
magnets, outside the ring but otherwise identical to the others. The current sup- 
plied to these two test magnets is in series with the beam magnets to assure the 
same current, and the magnetic field is measured with nuclear magnetic resonance 
(NMR) probes. The naive beam energy, 

-E-beam = electrou charge x magnetic field x CESR radius (3.1) 

is then corrected for shifts in the RF frequency, steering and focusing magnet 
currents, and the voltage of the electrostatic separators. This full reckoning misses 
the true beam energy by 0.172%, which is 18 MeV in -Ecm near the T masses, 
but it is very stable with time and tracks beam energy differences with the same 
precision. Such a beam energy measurement will not improve the world knowledge 
of T masses, but the scale uncertainty is small enough to have negligible impact 
on the width, and therefore the area, of the resonance scans. 

Distributions of electron and positron energies in the CESR beam are 0.057% 
wide at the T(IS') (this is the ratio of the standard deviation to the mean) and this 
width scales linearly with beam energy. The beam energy distribution is Gaussian. 
Our lineshape data, which is the world's most sensitive test of e^e~ beam energy 
distributions near 10 GeV, show no deviations from a pure Gaussian distribution. 

The beam energy spread can vary by as much as 1-2% a month, due to pertur- 
bations in the beams' orbits from environmental conditions. We observed such a 
shift (Figure 3.2), coincident with large corrections to the horizontal steering mag- 
nets to compensate for the new orbit. We use records of these changes to track 
potential shifts in the beam energy spread, and allow shifts in the widths of the 
lineshapes by adding fioating parameters to the fit. 

The beam-beam interaction region is a ribbon 0.18 mm tall (out of the ring 
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Figure 3.1: Greeting card from CESR in 1979, demonstrating its success in colliding 
e"'"e~ at 10 GeV by scanning the lineshapes of the T(IS'), T{2S), and T{3S) (left 
to right). 
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Figure 3.2: Beam energy spread fit results in the center-of-mass {AE = v2 single 
beam energy spread) divided by Eqm^, which is nearly constant, as expected. The 
third T(15') measurement (April 2002) is 2.4% lower than the second (March 2002). 
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plane), 0.34 mm wide (in the ring plane but perpendicular to the beam axis), and 
1.8 cm long (in the ring, along the beam axis). Electrons and positrons may collide 
anywhere within this constrained distribution, and its center drifts by about 4 mm 
along the beam-line and 1-2 mm perpendicular to it on a monthly timescale. We 
can easily track these changes. 

In addition to e+e^ collisions, beam particles can interact with gas nuclei inside 
the beam-pipe and with the wall of the beam-pipe itself (2.1 cm in radius). To 
minimize the number of beam-gas events, the beam-pipe is kept evacuated at 
2-4x10^^ torr. Beam-wall events are minimized by focusing. The electron and 
positron currents can vary independently, so electron-induced beam-gas and beam- 
wall and positron-induced beam-gas and beam-wall rates are not identical. In 
fact, we find that positron-induced rates are typically twice electron-induced rates, 
suggesting a difference in cross-sections. 

We collected data in eleven dedicated scans of the T(IS') and one high-energy 
point (100 MeV above the T(IS') mass), totalling 0.09 fb~^, and added to this 
18 fb^^ of subsequent on-resonance peak data (adjacent in time and limited to 
48 hours after the beginning of the scan). We obtained six T(2S') scans with 
a high-energy point (60 MeV above the T(2S') mass), totalling 0.05 fb^^ and 
added 0.03 fb^^ of subsequent peak data, and seven T(3S') scans with a high- 
energy point (45 MeV above the T(3S') mass), totalling 0.08 fb^^, and added 
14 fb~^ of subsequent peak data. We present all the individual scans in Table 3.1. 
In addition to scan data, we used the full 0.19 fb^^, 0.41 fb^^, and 0.14 fb^^ 
off-resonance datasets, 20 MeV below the T(IS'), T(2S'), and T(3S') masses, to 
subtract continuum backgrounds. 
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Table 3.1: Scans of the T(IS'), including associated on-resonance peak data. 
"Spread" indicates groups of scans which have the same beam energy spread (same 
labels as in Figure 3.2). The last entry is a point taken 100 MeV above the T(IS') 



mass. 



Scan 


Int. Lum. (pb ^) 


Run Range 


Dates 




Spread 


janl6 


6.7 


123164-123178 


Jan. 15-16, 


2002 


T(15)i 


jan30 


52.7 


123596-123645 


Jan. 30-Feb. 1, 


2002 


T(15)i 


feb06 


26.3 


123781-123836 


Feb. 6-8, 


2002 


T(15)i 


febl3 


7.8 


124080-124092 


Feb. 19-20, 


2002 


T(15)i 


feb20 


21.0 


124102-124159 


Feb. 20-22, 


2002 


T(15)i 


feb27 


23.9 


124279-124338 


Feb. 27-Mar. 1 


,2002 


r{is)2 


mar 06 


19.6 


124436-124495 


Mar. 6-8, 


2002 


TilSh 


mar 13 


25.9 


124625-124681 


Mar. 13-15, 


2002 


r{ish 


apr08 


7.2 


125254-125262 


Apr. 8-9, 


2002 


T{lSh 


apr09 


5.6 


125285-125295 


Apr. 9-10, 


2002 


T(15)3 


aprlO 


42.3 


125303-125358 


Apr. 10-12, 


2002 


T(15)3 


+ 100 MeV 


11.6 


124960-124973 


Mar. 27-28, 


2002 
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Table 3.2: Scans of the T{2S) and T(3S'), including associated on-resonance peak 
data. "Spread" indicates groups of scans which have the same beam energy spread 
(same labels as in Figure 3.2). The "+60 MeV" and "+45 MeV" entries are points 
taken 60 and 45 MeV above the T(15') and T(3S') masses, respectively. The 
dates of the "+45 MeV" data-taking overlap the "dec26" scan, but the integrated 
luminosity we quote do not. 



Scan 


Int. Lum. (pb^^) 


Run Range 


Dates 




Spread 


may29 


14.6 


126449-126508 


May. 29-3 L 


, 2002 


T{2S) 


junll 


9.9 


126776-126783 


Jun. 11-12, 


2002 


T{2S) 


junl2 


23.6 


126814-126871 


Jun. 12-14, 


2002 


T{2S) 


jullO 


18.8 


127588-127615 


Jul. 10-11, 


2002 


T{2S) 


jul24 


5.8 


127924-127933 


Jul. 23-24, 


2002 


T{2S) 


aug07 


9.3 


128303-128316 


Aug. 7-8, 


2002 


T{2S) 


+60 MeV 


4.9 


127206-127219 


Jun. 26-27, 


2002 




nov28 


27.5 


121884-121940 


Nov. 28-30, 


2001 


T{3S)i 


dec05 


41.3 


122069-122126 


Dec. 6-8, 


2001 


r{3S)2 


decl2 


41.3 


122245-122298 


Dec. 12-14, 


2001 


r{3S)s 


decl9 


24.2 


122409-122452 


Dec. 19-22, 


2001 


T{3S), 


dec26 


27.7 


122535-122579 


Dec. 25-26, 


2001 


r{3S), 


jan02 


27.7 


122766-122821 


Jan. 2-4, 


2002 


T(35)6 


jan09 


44.5 


122993-123044 


Jan. 9-11, 


2002 


T(3^)7 


+45 MeV 


10.8 


122568-122575 


Dec. 26, 


2001 
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3.2 CLEO Detector 

The CLEO detector is a general-purpose assembly of detectors built concentri- 
cally around the CESR interaction point [10] [11]. This analysis uses only three of 
CLEO's detectors: the silicon vertex detector and central drift chamber for iden- 
tifying charged particles, and the Csl crystal calorimeter for measuring electron 
and photon energies, and for simple particle identification. The CLEO-III appa- 
ratus, which is the generation of CLEO in operation in 2001-2002, is depicted in 
Figure 3.3. 

We define the z axis of our coordinate system to be parallel with the beam-line, 
pointing in the direction of the incident positron current (west). Our coordinate 
system is right-handed, with y pointing up and x pointing away from the center 
of the CESR ring (south). The origin of the coordinate system is at the center 
of the drift chamber, and lies within 1-2 mm of the beam-beam collision point. 
This coordinate system is illustrated in Figure 3.4. The CLEO detector has an 
approximate cylindrical symmetry around z, so we also define the polar angle 
^ of a particle trajectory originating at the origin to be the angle between the 
trajectory and the beam-line, or ^ = tan^^ ( A/aJM-y^Z-z ) . We often use cos^ and 
cot Q to describe the polar angle. The azimuthal angle is the angle for which 



COS0 = xj \/ x^ + ip and sin0 = yj^x'^ + y'^. 

The silicon vertex detector and the drift chamber both detect tracks left by 
charged particles by collecting charge left in the wake of ionizing, high-energy 
particles. In the vertex detector, the ionized medium is silicon, cut into strips 
held perpendicular to the trajectories of most particles (Figure 3.5). The charge 
is conducted out of the detector for amplification along traces which are parallel 
to the beam-line on one side of the strip and perpendicular to it on the other, so 
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Figure 3.3: An isometric view of the CLEO-III detector, highlighting the sihcon vertex detector, the drift chamber, and the 
Csl crystal calorimeter in gray. 
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that the two-dimensional point of intersection may be reconstructed. 

In the drift chamber, high-energy charged particles ionize a helium-propane gas 
(60% He, 40% CsHg) in a strong electric field generated by wires strung across the 
detector volume, parallel with the z axis. One quarter of these wires, called sense 
wires, are held at +2100 V, and the remaining three quarters, called field wires, 
are held at ground. The resulting field causes the freed electrons to drift away 
from the field wires toward the sense wire, which conducts the charge to amplifiers 
for analysis (Figure 3.6). As the electrons drift several millimeters toward the 
sense wire, they ionize more atoms, causing an avalanche that provides a 10^ 
amplification. We measure the time between the first ionization (estimated from 
bunch collision times) and charge collection on the sense wire to reconstruct the 
distance of closest approach of the high-energy charged particle to the sense wire, 
through the known electron drift speed of 28 /im/ns. This technique provides an 
average resolution of 88 /xm in the x-y plane. Sensitivity to z position is obtained by 
tilting the outer wires, presented in more detail in Figure 3.7. The outer 31 layers 
of wires, called the stereo section, are tilteded 21-28 mrad, yielding a z position 
resolution of 3-4 mm at each wire. The inner 16 layers, called the axial section, 
are untilted. 

Both tracking volumes are permeated by a 1.5 T magnetic field, pointing along 
the z axis. Charged particle trajectories are helical in this field: projections onto 
the x-y plane are circles. The polar angle Q of such a helical trajectory is a constant 
of the motion, but not 0. We measure the charge x momentum of particles through 
the radii of curvature of their tracks. Only electrons, muons, pions, kaons, protons, 
and deuterons are sufficiently stable and abundantly produced to be observed as 
tracks, and all of these particles have ±1 units of charge, so the radius of curvature 
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Figure 3.4: Coordinate system of the CLEO detector: (a) orientation with respect 
to the CESR ring and the Earth, (b) definition of 9 and (f). 



\ ^ 




Figure 3.5: Sihcon vertex detector geometry, projected onto the x-y plane. The 
trajectories of central [Q ~ 7r/2) collision products (dashed lines) are roughly 
perpendicular to the wafers of silicon. Triangles represent diamond rods holding 
the wafers in place. 
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Figure 3.6: Charge collection and multiplication in the drift chamber. 



a track pointing toward the FRONT of the detector 




tilted counter-clockwise 
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back ► FRONT 



Figure 3.7: Using tilted wires to obtain z information in the drift chamber. Tilted 
wires extrude lines in the x-y projection; position along a tilted wire indicates 
the z of the track helix near that wire. The closest wires to the track (in three 
dimensions) are highlighted. 
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provides access to monientum. The momentum resolution, dominated by drift 
chamber measurements, is 0.9% for beam-energy tracks, and position resolution 
near the interaction point, dominated by silicon vertex measurements, is 40 /xm in 
x-y and 90 /im in z. 

The outer radius of the drift chamber is 80 cm from the interaction point, and 
the outer edges are ±110 cm in z. The drift chamber's z range is more limited 
for smaller radii to accommodate the focusing quadrapole magnet, as shown in 
Figure 3.8. It will later be to useful to know that charged particles with more than 
60 MeV of z-momentum exit the detector before completing one half-orbit in the 
magnetic field. Such particles cannot generate multiple tracks by spiralling inside 
the detector volume. 

The Csl crystal calorimeter is sensitive to photons as well as charged particles, 
by presenting a transparent, high-Z material for them to interact with Electromag- 
netically. (Our thallium-doped Csl has a radiation length of 1.83 cm.) Incident 
electrons and photons are destroyed by this interaction, and replaced by a shower 
of equal total energy in less energetic photons, electrons, and positrons. Other 
particles deposit only a fraction of their energy. Visible light from the shower is 
collected on the back of the 30 cm-long crystals, from which the incident energy is 
reconstructed. Electrons and photons with energies near the 5 GeV beam energy 
are fully reconstructed with 1.5% resolution, but the energies of other particles is 
underestimated. Muons, for instance, deposit only 200 MeV in the calorimeter, 
regardless of incident energy. Combining the energy of calorimeter showers with 
track momenta is sufficient to identify and distinguish e"'"e~, /i'''/^", and 77 events 
with negligible backgrounds. 

The calorimeter geometry is composed of three parts: a barrel surrounding the 







Figure 3.8: A quarter of the CLEO detector, highlighting the sihcon vertex detector, the drift chamber, and the Csl crystal 
calorimeter in gray. Radial and z distances are in centimeters; the angular reach of the calorimeter barrel, the drift chamber, 
and the calorimeter endcap are illustrated by dashed lines. 
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tracking volume and two endcaps, beyond the tracking volume in z (see Figure 3.8). 
The calorimeter barrel covers polar angles with |cos^| < 0.85 and the endcaps 
extend this range to | cos^| < 0.97. The angular limits of the tracking volume is 
between these two: | cos^| < 0.93. 

When a threshold amount of activity is observed in the drift chamber and 
calorimeter, readout electronics are triggered to acquire a snapshot of the detec- 
tor and record all signals as an event. This activity is quantified in terms of the 
number of observed tracks and the number of showers above given energy thresh- 
olds. For speed in triggering, tracks are counted using a lookup table of drift 
chamber hits, trained by a simulation, and showers are approximated by summing 
calorimeter barrel output over 2x2 tiles, called clusters, and counting the number 
that surpass a given threshold. The number of AXIAL tracks is the number of tracks 
reconstructed in the axial section of the drift chamber, and STEREO is the number 
of tracks which can be extended into the stereo section. A CBLO cluster exceeds 
150 MeV, a CBMD exceeds 750 MeV, and a CBHI exceeds 1500 MeV. Real showers 
can be distributed over as many as four tiles, sometimes dividing their energy such 
that none of the clusters reach a threshold. This is a source of trigger inefficiency 
for final states that rely on shower information (Figure 3.9). After the data have 
been recorded, we reconstruct tracks and showers with much finer precision. 

We use several triggers to accept events, all of which are minimum-thresholds: 
an event is never rejected for having too many tracks or clusters. All of these 
triggers are active, and when an event is recorded, it is tagged with the names of 
the triggers it satisfied. The trigger relevant for this analysis are 

• two-track, which requires >2 AXIAL tracks, prescaled by a factor of 19 
(5.3% of the events satisfying this criterion are accepted). 
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Figure 3.9: The efficiencies of CBLO, CBMD, and CBHI identification as a function 
of fully-reconstructed shower energies for isolated showers. The efficiency curves 
are not symmetric around their thresholds because shower energy may be divided 
among tiles. 
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• hadron, which requires >3 AXIAL tracks and >1 CBLO, 

• rad-tau, which requires >2 STEREO tracks and (>1 CBMD or >2 CBLO), 

• e^-track, which requires >1 AXIAL track and >1 CBLO, and 

• barrel-bhabha, which requires 2 CBHI clusters on opposite sides of the 
calorimeter barrel. 

To count hadronic T decays, we select only those events which satisfied hadron, 
rad-tau, or e^-track, the three triggers that are efficient for hadronic decays. 
(This simplifies our efficiency study.) Note that a minimal condition for these 
three trigger is that at least one AXIAL track and one CBLO were observed. This 
minimal requirement is exact because STEREO tracks, being extensions of AXIAL 
tracks, are always less numerous than or equal in number to AXIAL tracks, and 
CBMD clusters are also CBLO clusters. 

Electron and positron beams are circulated in CESR for about an hour before 
their currents are exhausted from collisions. Data collected during this time is 
called a run, and is given a unique, ascending 6-digit identifier. Runs are the basic 
unit of CLEO data samples; in lineshape scans, we generally took one run at each 
Eqm point at a time. 

For some studies, we must simulate our entire detector on a computer. Such 
Monte Carlo simulations are most important in determining the efficiency-corrected 
cross-section for Bhabhas in CLEO, which is needed to measure the integrated lu- 
minosity of our datasets. While the total Bhabha cross-section is infinite, the 
efficiency-corrected cross-section, defined by observed Bhabhas, is finite and must 
be calculated theoretically. This calculation has two ingredients, the differential 
cross-section as a function of 6 and CLEO's efficiency for Bhabhas as a function of 
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9. The first ingredient is calculated with perturbative Quantum Electrodynamics 
(QED), but the second requires specific knowledge of our detector. While this 
efficiency may be approximated as a step function, in which CLEO observes all 
Bhabhas within a 9 range and misses all Bhabhas outside of this range, such a 
simplification would be bought at a high price in accuracy. For the 1% precision 
demanded by this analysis, we must consider all effects: detector geometry, elec- 
tron propagation and scattering in materials, sub-component response efficiency, 
fringe magnetic fields from the CESR magnets, et cetera. Our Monte Carlo simu- 
lation is based on the GEANT framework [14], and is carefully tuned to reproduce 
the real detector's output at all levels of analysis. 



Chapter 4 



Backgrounds and Event Selection 

To define a set of liadronic T decays, we will accept only those events which satisfy 
given criteria, or cuts. We want this set of events to include as many hadronic 
T decays as possible, to minimize the efficiency correction for lost events. We 
therefore only seek to reduce the backgrounds to a manageable level by cutting 
out regions of parameter space where the hadronic T contribution is minimal. We 
accomplish this with a set of four explicit cuts. 

With such an approach, we cannot completely eliminate backgrounds, espe- 
cially because continuum qq final states are identical to 8.9% of T decays. Instead, 
we estimate and subtract the backgrounds that remain after cuts, which we can 
do very accurately using control data. If we can accurately subtract any residual 
backgrounds after cuts, why cut at all? There are two reasons: large background 
subtractions introduce large statistical uncertainties, and the trigger itself selects 
events in a way which can be hard to predict, leading to systematic uncertainties. 
By imposing more restrictive event selections with fully-reconstructed data, we can 
render the trigger biases insignificant. 

4.1 Suppressing Backgrounds with Event Selection 

Bhabha scattering is our largest potential background before cuts, and among our 
largest backgrounds after cuts. We suppress Bhabhas by requiring the largest track 
momentum, |pmax|, to be less than 80% of -Ebcam- According to Monte Carlo, this 
rejects 0.15% of hadronic T events, but 99.73% of T ^ e+e^ and yU+/i~ (Fig- 
ure 4.1). We do not make a similar requirement on calorimeter shower energy. 
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which also peaks at -Ebeam for each electron in the Bhabha event, because we find 
track momentum measurements to be more stable in time than shower energy mea- 
surements. Unlike track momentum, which is a geometric measurement of wires 
in space, shower energy depends sensitively on the amplification of the calorime- 
ter read-out. This amplification is measured with 0.02-0.06% precision, but the 
Bhabha spectrum is so steep that 5% of Bhabha showers move across a reasonable 
threshold (75% of -Ebcam) with these fiuctuations in energy scale. 

It is also common to reject Bhabhas by requiring more than two tracks in the 
event: according to our simulation, 98.9% of hadronic events have more than two 
tracks and 99.4% of Bhabhas in the observable range have exactly two. We con- 
sidered this cut at an early stage in the analysis, but decided against it because we 
found the number of tracks distribution difficult to simulate for hadronic events. 
At that time, we intended to determine the cut efficiency with our Monte Carlo 
simulation, so this would have contributed significantly to the systematic uncer- 
tainty. Since then, we have found a way to measure hadronic efficiency without 
resorting to simulations, but we did not re-introduce the cut because we do not 
need it. The liVaxI cut reduces Bhabha contamination to approximately the same 
level as continuum gg, so the Bhabha contribution to the statistical uncertainty of 
the background-subtracted count is not dominant. (Several of our cuts imply that 
an event must generate at least one track, but this does not significantly affect the 
Bhabha background.) 

As previously mentioned, the cross-section of all continuum e^e~ -^ ff pro- 
cesses (/ is any fermion) fall off as 1/s while the two-photon fusion cross-section 
(e^e" -^ e^e~X) increases as log s. Continuum // may therefore be estimated and 
subtracted collectively, while two-photon fusion must be handled separately, pos- 
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Figure 4.1: Largest track momentum (|pmax|) in each event for data (solid his- 
togram) and Monte Carlo T decays (dashed), with all other cuts applied. Data 
have large backgrounds from Bhabhas and radiative Bhabhas. Hadronic T decays 
peak at 20% of -Ebcam, while T -^ e^e~ and ^i^ ii~ peak at 100% of -Ebcam- 
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sibly introducing systematic error if it is large. We therefore suppress two-photon 
fusion events by requiring the visible energy of the event, -Evis, to be greater than 
40% of Ecu- Visible energy is the sum of all track energies (determined from mo- 
mentum, assuming the charged particle to have a mass of 140 MeV) and neutral 
shower energies (neutral showers must be at least 7 cm from all tracks). If all 
particles in an event are detected, -Evis ~ -^cm- As seen in Figure 4.2, hadronic T 
events peak in _Evis at 80% Ecu, and there is a peak of non-T events at 15% Ecu- 
At least two-thirds of the events in this low--Evis peak are two-photon collisions, 
in which one incident electron has taken most of the center-of-mass energy, unde- 
tected, down the beam-pipe. We know this because two-thirds of events with less 
than 30% visible energy contain one low-momentum electron, whose charge and 
direction are correlated with the incident beams, and a highly anisotropic distri- 
bution of shower energy, presumably from the boosted hadron system X. Decays 
of T^T~ cover a broad spectrum of -Evis, due to energy lost in one or two neutrinos, 
extending but not peaking below our cut threshold. 

Rejecting low-i^^is events also protects our hadron count from uncertainties 
associated with trigger thresholds. In our simulations, only 0.07% of hadronic 
events with _Evis > 40% Eqm fail to trigger, so any fluctuations in the electronics 
will be on this level. The -Evis threshold, situated in the flat minimum between the 
two-photon peak and the signal peak, is minimally sensitive to fluctuations in the 
two-photon background and the signal efficiency. Since only 0.82% of simulated 
hadronic T decays fail this cut, any fluctuations in this efficiency will be well under 
a percent. 

In addition to the two-photon peak at an i^vis of 20% Ecu, there is an excess 
of events with _Evis just above 50% Eqm (Figure 4.2). These events are likely to 
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Figure 4.2: Total visible energy (-Evis) in each event for data (solid histogram) and 
Monte Carlo T decays (dashed) with all other cuts applied. At least two-thirds of 
the peak at 15% of Ecu in data are two-photon fusion events, and the peak above 
50% of Eqm are likely to be radiative Bhabhas, missing an electron. See text for 
a more complete discussion. 
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be radiative Bhabhas (e'^e^ -^ 7e^e^) in which one of the two electrons is lost. 
They contain neutral calorimeter energy and an energetic electron whose charge 
and direction is correlated with the incident beams, like the two-photon fusion 
events. However, visible energy in two-photon collisions is expected to be much 
less than 50% of -Ecm- 

The number of background events from a continuum process is proportional to 
the integrated luminosity, just like the number of signal T decays. Their contri- 
bution to the apparent cross-section will therefore be purely a function of Ecu- 
The same cannot be said for backgrounds that are not the product of beam-beam 
collisions. Beam-gas and beam-wall rates are a function of the individual beam 
currents, the gas pressure inside the beam-pipe (for beam-gas) and the extreme 
tails of the bunch shape (for beam- wall). Cosmic rays are abundant in our detec- 
tor, and the number of cosmic ray events is only a function of time. Integrated 
luminosity, integrated current, and time are approximately proportional (within a 
factor of two), so a continuum subtraction largely removes these effects, but not 
entirely. 

We suppress beam-gas, beam-wall, and cosmic ray events by requiring the event 
to originate near the beam-beam crossing point. To select events originating near 
this point in an x-y projection, we require at least one track to extrapolate within 
5 mm of the beam-line. We define dxy as the distance of closest approach of 
the closest track to the beam-line, and reject events with dxv > 5 mm. Tracks 
extrapolated from the tracking volume are corrected for momentum loss in the 
beam-pipe and silicon detector, and the location of the beam-beam crossing point 
is measured independently for each run, using the first 500 hadronic events. The 
dxY distribution (Figure 4.3) is much narrower than our 5 mm threshold: only 



53 

0.1% of beam-beam collision events fail this cut. This allows for ~1 mm errors in 
the beam-beam intersection measurement, which is far larger than expected. 

Our dxY cut is extremely effective at rejecting cosmic ray events. Cosmic rays 
rain uniformly into the detector, generating a uniform background to c^xy, which 
extends to 25 cm with our triggers. Only cosmic rays that pass within 5 mm of the 
beam-line survive. In principle, beam-wall events should also be eliminated, since 
they are generated in the beam-pipe, 2.1 cm from the beam-line. However, beam- 
wall events contain several tracks, any one of which may project into the accepted 
dxY region (Figure 4.4). By placing our requirement on the closest track, we bias 
this background to peak within our accepted region, diluting the effectiveness of 
the cut. 

Beam-gas and beam-wall events originate along the beam-line and beam-pipe, 
extending beyond beam-beam collisions in z. Placing a requirement on the closest 
track to the 2;-collision point would be ineffective for the same reason as above; 
beam-gas and beam-wall events both have many tracks, and the probability that 
one of these would project into the signal region (which is several centimeters 
wide) is not negligible. Instead, we reconstruct the z position of the event vertex 
using all tracks, and call this quantity dz- The CLEO event vertexing algorithm 
is not useful because it was designed for signal reconstruction and fails to fit too 
many beam-gas and beam-wall events. Instead, we developed a simple algorithm 
of our own. Tracking resolution is such that most of the tracks from a beam-beam 
collision intersect within 0.1 mm of a common origin in the x-y plane, and the 
number of intersections near this point grows rapidly with the number of primary 
tracks. Accidental track intersections far from this point grow more slowly. We can 
therefore determine the event vertex very accurately by averaging the z positions 
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Figure 4.3: Distance from the closest track projection to the beam-hne for data 
(sohd histogram) and Monte Carlo T decays (dashed) with all other cuts applied. 
The sign is related to the orientation of the track's curvature. The flat background 
in data is due to cosmic rays. 
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Figure 4.4: Though beam-wall events are centered at the beam-pipe, tracks may 
still project within 5 mm of the beam-line, thus passing the dxy cut. 

of track-track intersections. We define the z position of an x-y intersection to 
be halfway between the z positions of the two track helices, evaluated at the x-y 
intersection point. If the intersection is a true three-dimensional vertex, the tracks' 
z positions will be nearly equal. We weight these intersections with uncertainties 
propagated from the track uncertainties, the tracks' z separation, and the x-y 
distance to the beam-line added in quadrature, to prefer true intersections from 
the primary vertex. We plot this dz distribution in Figure 4.5, and cut very loosely 
at 7.5 cm, to allow for errors in the beam-beam intersection measurement. 

It is also possible to use track intersections to distinguish beam-wall events 
from beam-gas. The distance of the closest track-track intersection to the beam- 
line will be nearly zero for beam-gas events, but peak below the beam-pipe radius 
for beam-wall events (because selecting the closest intersection to the beam-line 
biases the distribution toward zero). In Figure 4.6, we plot the distribution of 
closest intersections for events with dxv < 5 mm from data with only one beam in 
CESR. We see that the dxv cut reduces beam-wall to the extent that it is approx- 
imately as common as beam-gas. A more sophisticated average of intersections 
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Figure 4.5: Location of the event vertex in z, according to our algorithm, for data 
(sohd histogram) and Monte Carlo T decays (dashed) with all other cuts applied. 
Data and Monte Carlo differ in the z-length of the beam-beam overlap region. The 
flat background in data is primarily beam-gas and beam-wall, though Monte Carlo 
indicates that 0.5% of T decays are misreconstructed and extend beyond the cut 
threshold. 
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could help to discriminate between beam-gas and beam-wall, but as we will see 
in Subsection 4.3.4, the two processes combined are a small contamination, about 
0.2% of the continuum for most runs. We therefore will not attempt to correct for 
beam-wall and beam-gas separately. 

4.2 Data Quality Requirements 

Not all data were collected under ideal conditions, so we applied some general 
criteria for rejecting bad runs. As the data were collected, two CLEO operators 
inspected the data for hardware failure. In the most serious cases, these data were 
eliminated from all CLEO analyses, but if the effect was limited, it was listed 
in a "bad runs" file (/home/dlk/Luminosity/badruns3S). We rejected any runs 
that were flagged with drift chamber, silicon vertex detector, or Csl calorimeter 
problems. 

We want a robust measurement of cross-section, and cross-section is constant 
with time, even as the beam currents are depleted during a run. We therefore 
checked for variations in cross-section during each run by comparing hadronic 
events and 77 events in hundredths of each run. This ratio fluctuates statistically, 
but we found two examples in which the drift chamber lost sensitivity to tracks 
before the calorimeter lost sensitivity to showers in the last few minutes of the run 
(Figure 4.7). Most likely, the drift chamber lost high voltage just before the end 
of the run. 

To catch more instances of this kind of failure, we also compared the rate 
of trackless Bhabhas to total Bhabhas. We recognize the e^e~ final state by 
the two beam-energy showers it produces in the calorimeter, curved 0.1 radians 
away from perfect collinearity by the magnetic field. Twenty-five runs had high 
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Figure 4.6: Closest track-track intersection to the beam-line in special data with 
one, non-colliding beam in CESR (all events are beam-gas, beam-wall, and cosmic 
rays). The rough peak below 0.6 cm is mostly beam-gas, and the broad peak from 
1 to 2 cm is due to beam-wall events. 
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Figure 4.7: Ratio of hadronic events to 77 events as a function of time through 
two runs. The hadronic cross-section drops to zero in the last 3% of these runs 
(the dashed hues are hnear fits to the data.) 
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trackless Bhabha rates (above 0.3%), and all of the trackless Bhabha excesses 
were in the same hundredth of a run (usually the last). Ten of these (presented in 
Figure 4.8) were crucial to the resonance scans and therefore not rejected. Instead, 
we determined the cross-section from the first 99% of these runs. The twenty-seven 
runs with drift chamber failures are listed in Table 4.1. 

The 77 final state, which we use for some diagnostic checks, is accepted only 
by the specialized barrel-bhabha trigger. We studied the efficiency of this trigger 
with Bhabha events and discovered eleven runs with very low efficiency, which we 
rejected, though these failures would only have affected our cross-checks. They are 
also listed in Table 4.1. 

We also tested the quality of the drift chamber and calorimeter output by 
counting unphysically high-energy tracks and showers. In good data, less than 1% 
of Bhabhas will generate a track or a shower with momentum or energy above 
120% -Ebcam- In a contiguous block of data on March 7, 2002, the fraction of high- 
momentum tracks abruptly increased to 3%. We see that the Bhabha peak for 
these runs has a high-energy tail (Figure 4.9), which suggests that the momen- 
tum in a fraction of tracks is overestimated. If this hypothesis applies to tracks 
with lower momenta, events may fail the |pmax| cut due to anomalous momentum 
measurements, changing the cut efficiency. We exclude these runs. On a separate 
occasion, December 16, 2001, the rate of high-energy showers abruptly rose to 3%. 
In this case, we observed that most of the unphysical showers occupy a regular 
block in the calorimeter barrel, indicating a read-out issue (Figure 4.10). Only 
our i?vis cut depends on shower energies, and in particular, only showers that can- 
not be associated with any track, so the influence of calorimeter malfunctions on 
our hadronic efficiency is limited. However, we will use showers in the calorime- 
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Figure 4.8: Bhabha events with no observed tracks as a function of time through 
ten runs. In all of these cases, there is an excess in the last 1% of the run. 
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Table 4.1: Runs rejected for hardware/calibration reasons. 



Drift chamber failed at the end of the run 


121476 


121748, 


121822, 


121847, 




122685 


123281, 


123411, 


123436, 




123847 


123873, 


124816, 


124860, 




124862 


125367, 


126273, 


126329, 




127280 








barrel-bhabha trigger inefficiency 


121928 


121929, 


121953, 


127951, 




127955 


130278, 


121710, 


121930, 




121944 


121954, 


123884 




Overestimated track momenta 


124452 


124454, 


124456, 


124458, 




124462 


124464, 


124465, 


124466, 




124467 


124469, 


124472, 


124473, 




124474 


124475, 


124477, 


124478, 




124479 


124480 






Overestimated barrel shower energies 


122331 


122335, 


122336, 


122339, 




122341 


122342, 


122344, 


122345, 




122349 


122350, 


122352 




Large cosmic ray /beam-gas backgrounds 


122353 


126341, 


129522 




Large, unidentified backgrounds 


121595 


122093, 


122330, 


126510 


Too little data for tests 


123013 


123014 
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ter barrel to identify Bhabhas and 77 events, so we exclude these runs as well. 
Another calorimeter malfunction, this time in the endcap, occured on December 
25-29, 2001. The -Evis spectrum for off-resonance runs in this time period is not 
distorted by excess background from a high-side tail on the two-photon fusion peak 
(see Figure 4.11), so we do not exclude these runs. All rejected runs are listed in 
Table 4.1. 

We rejected a handful of runs due to high background rates. From Figure 4.19, 
we set a 5% upper limit on acceptable cosmic ray yields relative to the continuum 
yield, and an upper limit of 2% on beam-gas. Three runs failed these criteria. We 
also noticed that the fractions of hadronic, Bhabha, 77, and /i+yU~ events dropped 
abruptly in the middles of four runs, indicating a sudden turn-on of some large 
background. We rejected these, too. Finally, two runs had so little data (16,695 
events total) that it was difficult to perform any of the above tests. We rejected 
them for convenience. 

This analysis combines small "scan" datasets, taken on the T resonances but 
not at its maximum, with off-resonance and "peak" data taken at the maximum 
cross-sections. The scan data were acquired specifically for this analysis and there- 
fore were not rejected lightly. (Only one run in Table 4.1 is a scan run: 124452.) 
The peak data are less valuable, and even after the selections described above, 
far more is available than is necessary. A measurement of the area of an T line- 
shape (i.e. Fee) can be conceptually decomposed into width measurements and 
height measurements, in which the fractional uncertainty in the area is the sum 
of the fractional uncertainty in the width and in the height, in quadrature. Scan 
data constrain both the width and the height, while peak data constrain only the 
height. Adding peak data to a fit will always reduce the statistical uncertainty. 
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Figure 4.9: The largest track momentum in off-resonance T(IS') (solid histogram) 
and March 7, 2002 runs (points with errorbars). 
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Figure 4.10: The locations of unphysical shower energies on December 16, 2001. 




Figure 4.11: The visible energy spectrum in off- resonance T{3S) (solid histogram) 
and December 25-29, 2001 runs (points with errorbars). 
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though this reaches an asymptotic hmit as the uncertainty comes to be dominated 
by the width measurements. However, as the beam energy cahbration drifts with 
time, cross-sections shghtly off the peak of the resonance may be represented as 
being exactly on-resonance, thereby biasing the height measurement. We accepted 
no more peak data than what is necessary to bring the statistical uncertainties 
within 5% of their limiting values. Since we are concerned with potential drifts 
with time, we re-expressed this limit as a time limit: we only include peak data in 
a lineshape fit if this data were taken less than 48 hours after the beginning of a 
scan. We imposed no limit on off-resonance data. 

We rejected a T(IS') scan, acquired on April 3, 2002. This scan is missing 
key cross-section measurements on the high-energy side of the peak (Figure 4.12), 
which makes it difficult to assess uncertainties in the beam energy and the beam 
energy spread. This scan does include cross-section measurements well above the T 
mass, and may have been the victim of miscommunicated beam energy requests. 
(Requests are made relative to the T mass, and single-beam energies used by 
CESR differ from our center-of-mass energies by a factor of two.) Its exclusion 
from the T(IS') fit affects the fit result by 0.12% with no appreciable difference in 
uncertainty. 

4.3 Subtracting Residual Backgrounds 

Backgrounds remaining after our cuts are summarized in Figure 4.13. We will 
discuss each of these, and their subtractions, in the subsections that follow. 
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Figure 4.12: The April 3, 2002 lineshape scan, overlaid by a fit to all other T(IS') 
scans. No data significantly constrain the high-energy side of the peak. 

4.3.1 Backgrounds that Vary Slowly with Beam Energy 

After our cuts, radiative Bhabhas and continuum qq dominate the background, 
adding a fiat, 8 nb plateau below our three T peaks (18 nb, 7 nb, and 4 nb, 
respectively) in apparent cross-section versus -Ecm- All continuum processes except 
for two-photon fusion evolve as 1/s, so we include such a function in our lineshape 
fits. The magnitude of this term is determined independently for the T(IS'), T(2S'), 
and T(3S') by the large off-resonance samples taken only 20 MeV below each T 
mass. The 1/s curve is the dashed line near the top of Figure 4.13. 

The first correction to the background curve is to add lower-energy T reso- 
nances, which have a l/(^/s — Mx) distribution. The magnitude of an ISR tail 
is set by the magnitude of the T resonance. We therefore fit T(IS'), T(2S'), and 
T(3S') in ascending order to obtain tail corrections from the previous fits. The 
T(IS') and T(2S') ISR tails under the T(3S') peak are labeled in Figure 4.13, and 
Figure 4.14 shows T(2S') and T(3S') off- resonance cross-sections with and without 
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Figure 4.13: The T{3S) lineshape in log scale, to illustrate backgrounds. The top 
dashed curve represents the sum of all backgrounds, which is dominated by 1/s 
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corrections included in "all backgrounds." Dashed curves represent ISR tails from 
charmonium resonances which are included in the two-photon fusion curve. The 
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this tail correction. 

To parameterize the logs correction for residual two-photon fusion, we fit the 
three off-resonance cross-sections to A/ s+B log s and present this fit in Figure 4.15. 
We find (8.0 ± 0.5)% of the apparent cross-section at 9 GeV to be due to the 
logs component. To see if this is plausible, we roughly estimate the two-photon 
background surviving our cuts by extrapolating the two-photon peak above our 
cut threshold in E^is (see Figure 4.15), yielding a two-photon fraction of 6%. This 
is consistent with our A/s + B\ogs fit. Other effects may contribute to part 
of the logs term, such as Eqm dependence in our cut efficiency for continuum 
events, a slow variation in the hadronic continuum cross-section, and ISR tails 
from charmonium resonances {J /ip and ip' , see Figure 4.13 for scale). All of these 
effects vary slowly with Equi so our parameterization for large differences in Eqm 
(900 MeV from T(IS') to T(3S')) applies to small differences in i^cM as we project 
the T(IS'), T(2S'), and T(3S') off-resonance cross-sections below each peak. The 
difference in cross-section between a pure 1/s curve and the fully parameterized 
curve is only 0.04% at the peak. 

4.3.2 Continuum-Resonance Interference 

As discussed in Section 2.1, resonant T ^ qq interferes with continuum qq. We 
must therefore also add a crint(-E'cM) term to our fit function (see Equation 2.8). 

In this analysis, we assume that e^e^ ^ qq ^ hadrons interferes with e^e^ -^ 
T — i> gg — i> hadrons but not e^e^ ^ T ^ ggg -^ hadrons, though the latter may 
share some final states which are indistinguishable from qq decays. (Interference 
from gg'-f is negligible because its branching fraction is only 3% of ggg and most 
gg'j events have a distinctive, high-energy photon.) For interference between qq 
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Figure 4.14: Off-resonance cross-section measurements versus Eqm, with and with- 
out ISR tail corrections. (Corrected data are at the centers of the circles.) The 
solid curve is the best fit to A/s + -Blogs, and the dashed curve is 1/s only, 
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70 



O 

d 

c 
> 

p-1 



6000 



5000 



4000 



3000 



2000 



1000 














0.2 



0.4 



0.6 



0.8 



1.2 



'VIS 



EcM for off-resonance data 



Figure 4.15: The visible energy (-Evis) of off-resonance data, overlaid with a fit to 
the low-i^vis peak;. The fit is Gaussian on the low-energy side and Lorentzian on 
the high-energy side, and is used to roughly estimate the two-photon fusion events 
which survive after the 40% of Eqm cut (dashed vertical line). 



71 

and ggg decays, quantum states must remain coherent through the hadronization 
process. This effect has been observed in J /ip and ip' -^ tt+tt^ and K^K^, but it 
is unclear if the effect is significant when summed over all final state amplitudes, 
since they may cancel. The phase difference between qq and ggg for the inclusive 
process is also unknown, and some phase differences cannot be constrained by our 
lineshape fits. We will therefore only assume parton-level interference, and discuss 
full hadronic interference as a fitting issue in Chapter 8. 

4.3.3 Backgrounds from T 

Since we are selecting hadronic T events, T -^ e+e^, fi^fi^, and t^t^ are back- 
grounds which peak under the hadronic T signal. We have no control sample for 
leptonic T modes, so we estimate these with a Monte Carlo simulation: negligible 
e+e^ and fi^fi^ survive the |pmax| cut (0.22% and 0.25%), even with final-state 
radiation (T -^ '~fe~^e~ and '~ffi~^fi~) modeled by PHOTOS. Our cuts and trigger 
are 57% efficient for t~^t~ , however. A tau lepton may decay into several hadrons, 
making it difficult to distinguish from hadronic T decays. Tau-pairs are rejected 
primarily by the E^is cut, as their visible energy spectrum is very broad due to 
neutrinos in the final state. 

We will need to subtract t^t~ events from the hadronic T count. The i?cM 
dependence of T ^ t^t^ is the same as T ^ hadronic, though the magnitudes of 
the resonant and interference terms both differ. The resonant t^t^ contribution 
is a factor of i3^^/i3had times smaller than the hadonic resonance, and the r^r^ 
interference term has a ctint of 0.20, 0.37, and 0.27 for the T(IS'), T(2S'), and 
T{3S), respectively. Continuum t~^t~, like continuum qq, is included in the 1/s 
term. When we estimate systematic uncertainties in the lineshape parameteriza- 
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tion in Section 8.2 (page 123), we will note that the uncertainty in Brr overwhelms 
the uncertainty in t^t~ efficiency and t^t^ ajnt, so only the branching fraction 
uncertainty must be propagated. 

4.3.4 Beam-Gas, Beam- Wall, and Cosmic Rays 

The non-beam-beam backgrounds are not a strict function of integrated luminosity, 
so we will need to explicitly subtract them from the hadronic T count for each run. 
To do this, we identify cosmic ray events, beam-gas, and beam-wall events in every 
run with special cuts. We then use control samples containing only cosmic rays 
or cosmic rays, beam-gas, and beam-wall to determine how to relate the number 
of non-beam-beam backgrounds that we counted to the number that survive our 
hadronic cuts. We then subtract this excess. 

To identify cosmic rays, we require the following. 

• No track may project within 5 mm of the beamspot (I^xyI > 5 mm). 

• The event must contain at least two tracks, since our track reconstruction 
algorithm identifies the descending-radius part of the cosmic ray as one track 
and the ascending-radius part as another. 

• The normalized dot product of the two largest track momenta {pi-p2/\pi\\p2\) 
must be less than -0.999 or greater than 0.999, since the angles of these two 
tracks differ only due to tracking resolution (though the orientation may be 
confused by hits with unexpected drift times). 

• The total calorimeter energy must be less than 2 GeV, consistent with two 
minimally-ionizing muon showers, and 
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Table 4.2: Run numbers for beam-gas, beam-wall, and cosmic ray control datasets. 

no-beam 128706 128736 128741 128748 

electron single-beam 126828 126920 126922 

positron single-beam 126785 

• E^is > 4% of EcM for less sensitivity to trigger thresholds. 

These cuts are, by design, much more efficient for cosmic rays than our hadronic 
cuts, but the number of identified cosmic rays and the number of cosmic rays 
contaminating our hadron count is proportional. To determine this constant of 
proportionality, we apply both sets of event selection criteria to a data sample 
acquired with no beams in CESR. (The no-beam runs are listed in Table 4.2.) 
Figure 4.16 superimposes cosmic ray candidates from this no-beam sample on cos- 
mic ray candidates from a large beam-beam sample, indicating a clear separation 
between cosmic rays and beam-beam collisions in dxv- We assume that all events 
in the no-beam dataset which pass our hadronic cuts are cosmic rays, so the desired 
constant is just a ratio of the cosmic ray count to the hadronic event count in this 
sample. The effective cross-section of cosmic rays are plotted with uncertainties in 
Figure 4.13. 

Beam-gas and beam-wall events are hard to distinguish from one another, but 
they are both small backgrounds which depend on the electron and positron beam 
currents. This dependence is not identical, since beam-gas rates are proportional 
to the gas pressure inside the beam-pipe while beam- wall is not. However, the 
contamination from beam-gas and beam-wall combined is typically 0.2% of the 
continuum. Furthermore, our beam-gas and beam-wall cuts have a small back- 
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Figure 4.16: Distance from the beam-line (|(ixy| in logx-scale, where 1 = one 
meter) with all other cosmic ray cuts applied. Data with beam-beam collisions 
are the solid histogram with a peak at 0.0001 (0.1 mm) from collisions and 0.01 
(10 cm) from cosmic rays. Data from the no-beam sample are the points with error 
bars, normalized to equal numbers of cosmic rays. The dashed vertical line is the 
cut boundary at 5 mm. Triggers fail to accept cosmic rays beyond 25 cm. 
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ground from beam-beam data, meaning that our estimate is too large. Instead of 
subtracting all of this estimate, we inflate our uncertainty. 

To identify beam-gas and beam- wall events (which we will call beam-nucleus), 
we require 

• Mxy| < 5 mm, \dz\ > 7.5 cm, 

• Ipi ' 'P2\/\Pi\\P2\ < 0.9 to further reject cosmic rays, 

• at least two tracks, and E^is > 4% of -Ecm- 

To distinguish between electron-induced beam-nucleus and positron-induced beam- 
nucleus events, we also cut on the net z-momentum of all tracks (p^z)- Fo^ ^ 
positron-induced event, we require p^^ > 10% of -Ebcam because incident positron 
momentum is in the positive z direction (see Figure 3.4). Electron-induced events 
must have p^J < —10% of -Ebeam- 

To relate the number of identified beam-nucleus events to the number that 
contaminate our hadronic event count, we employ data samples acquired with 
only one beam in CESR (also listed in Table 4.2). To use these samples, we must 
first subtract the cosmic rays using the technique described above. Figure 4.17 
demonstrates the separation of electron- and positron-induced beam-nucleus by 
their net z-momenta. This Figure also indicates that a small fraction, perhaps 10%, 
of our beam-nucleus candidates are contaminated by beam-beam events, probably 
two-photon fusion with a misreconstructed dz- The potential for contamination 
is also evident in Figure 4.18. Therefore, a beam-nucleus correction in analogy 
with the cosmic ray correction would be an over-subtraction of about 10%. The 
beam-nucleus estimates are typically only 0.1% of the continuum (Figure 4.19), 
so we subtract 50% ± 50% of the electron- and positron-induced beam-nucleus 
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estimates. The effective cross-section of beam-nucleus estimates are plotted near 
the bottom of Figure 4.13. The cosmic ray and beam-nucleus estimates for every 
run we used to determine Tee are plotted in Figure 4.19. 
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Figure 4.17: Net z-momentum of all tracks (p^^) with all other beam-nucleus cuts 
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Figure 4.18: Distance of the z event vertex from the center of the beam-beam 
distribution {\dz\ in loga;-scale, where 1 = one meter). Data with beam-beam 
colhsions are the sohd histogram with a peak above the plot window at 0.01 (1 cm) 
from beam-beam colhsions and a peak at 0.1 (10 cm) from beam-nucleus collisions. 
Data from the single-beam samples are the points with erorr bars, normalized to 
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Figure 4.19: Total counts of non-beam-beam backgrounds as a fraction of the 
continuum level and a function of run number. Dashed vertical lines separate 
T{3S), T(IS'), and T(2S') data-taking periods. Note that all three plots have 
different vertical scales: cosmic rays (a) are the most abundant, and the proportion 
of positron-induced beam-nucleus events (c) are typically twice that of electron- 
induced events (b). 



Chapter 5 
Hadronic Efficiency 

5.1 Motivation for the Data-Based Approach 

Inefficiency is in some sense the opposite of the problem of backgrounds: after 
removing the events which should not be in our hadronic T sample, we need to add 
in the events that are missing. In this Chapter, we will determine the probability 
that a hadronic T decay is included in our count, for each of the three resonances. 
These efficiencies are high, about 97% for each resonance. 

Often, efficiencies are determined from Monte Carlo simulations. One simulates 
all known decay modes, and constructs an aggregate efficiency 

e = J]e,i3„ (5.1) 

i 

where e^ is the efficiency of each mode. We don't directly use this method for two 
reasons. 

a. Hadronic decays are the result of the hadronization of bare quarks and gluons. 
This is a non-perturbative process which is only empirically approximated 
by LUND/ JetSet in the Monte Carlo. If we assume a non-perturbative QCD 
model to determine a non-perturbative QCD parameter, we would introduce 
a circular dependence that would have to be quantified. 

b. Our definition of hadronic T decays includes potentially unknown modes 
whose efficiencies may be very different from the hadronic modes we simulate. 
For instance, it is possible that T decays into invisible wiMPs with zero 
efficiency or that unknown QCD resonances may enhance decays to Kl or 
neutrons, which fail our E^is cut with greater probability. 
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Instead, we take advantage of our 1.3 fb^^ sample of T(2S') decays to study 
T{2S) -^ 7r+7r^T(lS') transitions. The T(IS') mesons in these decays are produced 
nearly at rest, and decay as they would from direct e"'"e~ -^ T(IS'). However, the 
2S ^ IS cascade events additionally include two charged pions which may satisfy a 
trigger and cause the event to be recorded, regardless of how the T(IS') decays. As 
an extreme example, we can use this technique to collect events featuring invisible 
T(IS') -^ vv or WIMP WIMP decays, which would be impossible with a direct 
T(IS') sample. 

We exploit this broad access to T(IS') decays to measure the T(IS') efficiency. 
From T(2S') -^ 7r+7r^T(lS') cascades, we select a subset in which the tt+tt^ by 
itself guarantees that the trigger will accept the event, so that we know that the 
trigger did not rely on decay products of the T(IS'). Assuming that there is no 
correlation between the kinematics of the vr'^vr^ and the branching fractions of 
the T(IS'), this subset of tt^tt^ is accompanied by a generic set of T(IS') decays: 
T(IS') decay modes are represented in the data sample with the same proportions 
as in nature. We may then apply our cuts to the T(IS') decay products in our 
sample to determine the fraction which succeed. This is the efficiency. 

5.2 Hadronic Efficiency of the T(15) 

We have two goals for our event selection in this study, to identify T(2S') -^ 
7r^7r^T(lS') candidates and to choose tt+tt^ candidates which are sufficient to sat- 
isfy the trigger. The tt^tt^ in T(25') -^ 7r+7r^T(lS') are kinematically constrained 
by the mass difference between the T(2S') and the T(IS'), so the mass of the system 
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Figure 5.1: Distinguishing the recoil mass of tt+tt" in kinematically-constrained 
T{2S) -^ 7i^7i^T{lS) events from other, accidental track-track combinations 
(backgrounds). The background distribution is not simple, but it has no struc- 
ture in the narrow region of our interest. 



recoiling against the two pions. 



m^ 



M- 



T(25) 



Vl^F+^«^ - Vl^F+^«^ ) -\Pl+P2 



(5.2) 



peaks at the T(IS') mass. This allows for excellent background rejection, because 
the peak from T(2S') -^ 7i^7i^T{lS) has a 3 MeV resolution while the background 
spectrum is much broader (see Figure 5.1). We require tt^tt" candidates to have 
a recoil mass (m^jr-rec) between 9.441 and 9.480 GeV. We further suppress back- 
ground by requiring the track helices that we identify as tt+tt^ to intersect in the 
x-y plane within 5 mm of the nominal beam-beam collision point. The track he- 
lices at this x-y point must also be less than 2.5 cm from each other in z, and their 
average z must be within 5 cm of the beam-beam collision point. The momenta 
used in Equation 5.2 are evaluated at the intersection point. 

To satisfy the trigger, we select tt+tt" track candidates with more than 150 MeV 
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of momentum perpendicular to the z axis (p±), so that their trajectories reach 
beyond the sixteenth layer in the drift chamber and satisfy the geometric require- 
ments for AXIAL tracks. A study of CLEO's trigger response to hadronic tracks 
revealed that the probability for a track with p_|_ = 150 MeV to be detected as an 
AXIAL track is 99.96%, and this probability grows with p_|_ [15]. We can therefore 
be at least 99.92% certain that events accepted by the two-track trigger (which 
requires two AXIAL tracks) with this cut on the vr^vr" did not rely on T(IS') decay 
products to be accepted. 

To simplify the process of excluding the tt+tt^ candidates when we apply our 
hadronic cuts, we additionally require each pion track to have more than 60 MeV of 
2;-momentum magnitude. As discussed in Section 3.2 on page 41, such trajectories 
exit the detector before completing one half-orbit in the magnetic field. This 
protects our sample from events with pions which spiral in the tracking volume, 
potentially generating many tracks, only one of which is identified as a pion from 
T(2S') -^ 7r+7r^T(lS'). With this cut, we can assume that each pion is responsible 
for only one track and the calorimeter showers associated with that track. 

There is occasionally more than one pair of tracks which satisfy these vr+vr" 
criteria in a single event. In this case, we choose a tt+tt^ candidate randomly. If we 
were to choose the tt+tt^ candidate that best reconstructs the T(IS') mass, would 
bias the non-cascade background to peak under the T(2S') -^ 7r+7r^T(lS') signal, 
complicating the background subtraction. Since the tt+tt^ candidate is chosen 
randomly, the background distribution has no structure on the scale of tens of 
MeV, and we can approximate it with a low-order polynomial. 

Selecting T(2S') -^ 7r+7r^T(lS') candidates in this manner, from events ac- 
cepted by the two-track trigger, and applying our hadronic cuts to the T(IS') 
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events in the peak, we find a hadronic T(IS') efficiency that is consistent with 
100% with 3% uncertainty, which is not satisfactory for our analysis. Hadronic 
efficiency is a factor in Feerhad/rtot, so the fractional uncertainty in hadronic effi- 
ciency adds to the reeFhad/rtot fractional uncertainty in quadrature. 

The main culprit is the prescaled two-track trigger, which randomly rejects 
94.7% of the events that satisfy its two AXIAL track criteria. We can circumvent this 
loss of data and improve the precision of our result by splitting the T(IS') hadronic 
efficiency into two factors. Define an T(IS') decay as "visible" if the trigger records 
one AXIAL track from the T(IS') decay products and one CBLO cluster from either 
the T(IS') decay or the tt+tt^. Define evis to be the probability that an T(IS') decay 
is visible, and ecuts to be the probability that a visible T(IS') decay is selected by 
our triggers and cuts. The hadronic T(IS') efficiency is the product of Cvis and 
Ccuts- We will determine evis, a number which is very close to 100%, using the 
two-track trigger, and ecuts using the hadron trigger, which selects visible T(IS') 
decays accompanied by tt+tt^. Since (1 — evis) is such a small inefficiency, a large 
fractional uncertainty in (1 — evis) propagates to a small fractional uncertainty in 

Our definition of "visible" is tailor-made for the hadron trigger in T{2S) -^ 
7r+7r^T(lS') events. The hadron trigger requires three AXIAL tracks and one CBLO 
cluster. We have selected 7r"''7r~ kinematics such that each pion must generate one 
AXIAL track, so the T(IS') decay is responsible for the third track and possibly 
the CBLO. Whatever the probability is that the CBLO comes from the tt+tt^, rather 
than the T(IS'), it is the same while measuring evis as it is while measuring ecuts- 
This gives us the freedom to determine the cut efficiency with a large set of events 
from the hadron trigger, knowing that we can correct for the bias it introduces 
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with evis- 

In our article in Physical Review Letters [16], we discuss this technique more 
succinctly using two factors, ehtrig and ecuts- In that article, ehtrig = Cvis, and ecuts 
has the same meaning in both. 

5.2.1 Determination of evis 

To determine evis, we apply the above procedure, replacing our full set of cuts 
for the "visible" condition we have just defined. The recoil mass distribution of 
TT+TT^ candidates accepted by the two-track trigger is shown in Figure 5.2(a). 
The real T(2S') -^ 7r+7r^T(lS') events peak near the T(IS') mass of 9.46 GeV, 
and backgrounds are distributed linearly underneath. To measure the fraction of 
T(IS') decays that are "invisible," we select from this sample events that fail the 
hadron trigger. These are plotted in Figure 5.2(b), revealing no significant peak. 
This indicates that nearly all of the T(IS') decays are visible. 

To quantify this statement, we fit the data in Figures 5.2(a) and 5.2(b) to the 
same function and extract the ratio of T(2S') -^ 7i^7i^T{lS) yields. We fit the full 
two-track dataset (Figure 5.2(a)) to a double Gaussian with a 1:4 ratio of areas 
and a linear background. When we fit the invisible-T(lS') data (Figure 5.2(b)), 
we fix the mean and widths, imposing the assumption that momentum resolution 
of the 7r+7r~ is independent of the visibility of the T(15'). This very reasonable 
assumption purchases much of the statistical precision in the measurement. The 
ratio of fit yields, which is the probability that an T(15') decay will be invisible, is 
(0.67 ± 0.62)%. 

The invisible T ^ z/z/ and T -^ wimpwimp we discussed earlier are in- 
cluded in (1 — Cvis), though the 0.05% T — > z/z/ branching fraction is much smaller 
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Figure 5.2: Recoil mass of tt^tt^ candidates for (a) all events accepted by the 
two-track trigger, and (b) only those in which the rest of the event is "invisible." 
True T(2S') -^ 7r+7r^T(lS') events are contained in the peak near 9.46 GeV, which 
is highly suppressed in (b), indicating that very few T(IS') decays are invisible. 
The solid curve is a fit to the data, and the dashed curve is the linear background 
contribution. 
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than the statistical uncertainty. This measurement does, however, place a model- 
independent upper bound on unknown, neutral decays of T(IS'), which may be in- 
teresting for constraining models that feature new neutral particles such as wiMPs 
or large branching fractions for all-neutral hadronic events. This is a constraint on 
neutral decays because our definition of "visible" requires at least one track with 
I cos^l < 0.93 and p± > 150 MeV. The T(IS') branching fraction to events with 
no such tracks is less than 1.01% at 90% confidence level. 

This ratio of fit yields represents the fraction of T(IS') decays that are invisible, 
including leptonic T(IS') decays. Leptonic decays, particularly e+e~ and fi^n~, 
are more likely to be invisible than hadronic decays, since their final state consists 
of only two particles which are geometrically back-to-back: if one lepton disappears 
down the beam-pipe, the other probably will do so on the other side. We determine 
the probability for leptons to be invisible from leptonic Monte Carlo simulations 
(which include the (1 -|- cos^) angular distribution for leptonic decays through a 
virtual photon), and this probability is (10.91 ± 0.01)%. We use Equation 5.1 to 
determine the hadronic visibility efficiency, Cvis, from the total visibility efficiency 
(el°l = 1 - 0.0067) and the leptonic visibility efficiency [e^l"^ = 1 - 0.1091). 



(1 - 3B„ 



(100.02 ±0.62)%. (5.3) 



Slightly more than half of this probability distribution is above 100%, which is im- 



possible for a real efficiency, so we truncate the part above 100% and normalize the 

-0.45 J /O. 



remaining distribution to obtain an asymmetric uncertainty: evis = (99.59 '^^■'^^^ 



5.2.2 Determination of e 



cuts 



We select tt+tt candidates in the same way to determine ecuts, except that we 
choose events accepted by the hadron trigger rather than the two-track trigger. 



Figure 5.3(a) presents the recoil mass of these 7i~^tt^ candidates, recoihng against 
visible T(IS') decays. We apply our hadronic cuts (but not the trigger require- 
ments) to this sample, excluding tracks and showers associated with the 7i~^7i~, 
and plot the recoil mass for these in Figure 5.3(b). These "cut-failure" events have 
a prominent peak at the T(IS') mass due to all the visible T(IS') events which 
failed our cuts. (Most of them are leptonic decays.) 

To obtain the ratio of cut-failure T(IS') events to all visible T(IS') events, we 
avoid the fit procedure because of the potential for systematic uncertainties in the 
fit parameterization. In the evis study, the statistical uncertainty in the number 
of invisible events was almost as large as the number of invisible events itself, and 
this overwhelmed any bias introduced by the fit function shape. Here, the number 
of cut-failures is significantly greater than zero and ought to be measured more 
precisely. Instead of fitting, we count events with a recoil mass between 9.454 and 
9.472 GeV and subtract the backgrounds, which have been determined by a linear 
fit to the sidebands (between 9.441 and 9.480 GeV, excluding the signal region). To 
estimate the systematic uncertainty in yield due to assuming a linear background, 
we repeat the procedure with the largest quadratic term allowed by the data. 

The ratio of cut-failure T(IS') events to all visible T(IS') events is (92.58 ± 0.13 
± 0.02)%, where the first uncertainty is statistical and the second is due to the 
parameterization of the background distribution. The leptonic modes account for 
most of this 7% inefficiency: correcting for leptonic modes as above (Equation 5.3) 
yields ecuts = (98.32 ± 0.21)%. 

The apparent T(IS') efficiency, measured in T(2S') -^ 7r+7r"T(lS') cascades, 
may differ from the efficiency of direct e~^e~ -^ T(IS') events because of the slight 
relativistic boost of the T(15') (7 = 1.005) or the potential for tt+tt^ showers to 
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Figure 5.3: Recoil mass of tt+tt^ candidates for (a) all events accepted by the 
hadron trigger, and (b) only those in which the rest of the event fails our cuts 
("cut-failure" events). The dotted vertical lines identify the signal region, and the 
dashed curve is a fit to the background, used to subtract background events from 
the signal region. 
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overlap T(IS') showers in cascade decays. To address these possibihties, we gener- 
ate hadronic Monte Carlo for the direct and the cascade cases, applying the same 
procedure to determine the hadronic efficiency from the cascade simulation. We 
observe no significant difference: the ratio of direct efficiency to cascade efficiency 
is 1.0014 ± 0.0022. We apply this as a correction primarily to propagate the 
uncertainty in this study. 

Incidentally, the cascade Monte Carlo prediction of Ccuts, (98.54 ± 0.22)%, 
agrees with our data-based measurement of (98.32 ± 0.21)%. We would have been 
correct, if not justified, if we had derived our cut efficiency directly from the Monte 
Carlo. We also extract |pmax|) -^visj dxYi and rfz from our cascade data and cascade 
Monte Carlo, and find that they agree fairly well (Figure 5.4). 

One correction is still missing from the ecuts we have derived: ecuts niust include 
the trigger efficiency, but we excluded trigger requirements from our measurement. 
While we can remove the tt^tt^ tracks and showers from our fully reconstructed 
data, it would be impractical to apply an analogous procedure on our trigger data 
for technical reasons. We therefore use Monte Carlo to determine the efficiency of 
the trigger once our cuts have been applied, a value of 99.87%. Figure 5.5 overlays 
data and Monte Carlo AXIAL, STEREO, CBLO, and CBMD distributions in direct T(IS') 
decays, showing fairly good agreement. The predicted inefficiency of 0.13% can 
therefore be trusted within 100% of itself, so we conservatively assign this as the 
systematic uncertainty. 

After all of these corrections, the efficiency of visible hadronic T(IS') decays is 
ecuts = (0.9832)(0.9987)(1.0014) = (98.33 ± 0.33)%. The efficiency of all hadronic 



T(IS') decays is evis x Ccuts = (97.93 



-0.56 J /O. 
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Figure 5.4: Our four cut variables, as seen in background-subtracted T(2S') -^ 
7i^7i^T{lS) events. Points with errorbars are data, the sohd histograms are 
T(2S') -^ 7i^7i^T{lS) Monte Carlo simulations with the same procedure applied, 
and the dotted vertical lines are the cut thresholds, (a) Largest track momentum 
(IPmaxI), divided by Mx(i5)/2, the equivalent of -Ebeam if this were a direct decay. 
The peak at 1 is due to e+e^ and /i+/i^ decays, cross-hatched in the Monte Carlo. 
(b) Visible energy (i^vis); divided by My(is), the equivalent of Eqm- (c) The dis- 
tance of the closest track to the beam-line (c/xy)- (d) The z- vertex of the event 
{dz). 
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Figure 5.5: The four cut variables used in hadron, rad-tau, and e^-track trigger 
decisions (defined on page 45). Points with errobars are continuum-subtracted 
data, histograms are Monte Carlo simulations, and the dashed vertical lines are 
cut thresholds. Trigger selections have been applied to data and Monte Carlo, 
though Monte Carlo without trigger selections are overlaid as dotted histograms 
(barely visible). 
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5.3 Hadronic Efficiency of the Y(25) and T(35) 

The T(3S') -^ 7r+7r^T(2S') and T{AS) -^ 7r+7r^T(3S') rates are too low to accumu- 
late large samples to study T(2S') and T(3S') efficiencies in the same way that we 
did T(IS'). Instead, we derive correction factors from the Monte Carlo that allow 
the T(IS') efficiency to be applied to the T(2S') and T{3S). 

The decays of the T(2S) and T{3S) differ from those of the T(IS') in two ways. 
The T(2S') and T(3S') decay products are slightly more energetic, as they originate 
in a state of higher Eqm, and the T{2S) and T{3S) can decay via transitions to 
lower bb resonances. The ffist correction is very small because our |pmax| and E^is 
cut thesholds are constant fractions of Eqm, and the difference in Eqm from T(IS') 
to T(35') is only 10%. In the Monte Carlo simulation, the T(2S') and T(35') 
efficiency for ggg, gg'~f, and qq is only 0.2% lower than the T(IS') efficiency. 

Most T(2S') and T(3S') transition decays have the same efficiency as ggg, gg'~f, 
and qq, so they have no impact on total hadronic efficiency. The exceptions are 
transitions that result in a lower T resonance decaying into e~^e~ or fi^fi". Accord- 
ing to the Monte Carlo simulation, these "cascade-to-leptons" decays have (0.69 ± 
0.22)% efficiency for T{2S) and (0.38 ± 0.19)% efficiency for the T(35)— almost 
zero. Therefore the efficiency correction will be approximately (1 — Bcas), where 
Bcas is the branching fraction for these modes. 

We determine this branching fraction from the data by counting cascade-to- 
leptons events relative to T ^ fifi in the full T(2S') and T(3S') datasets. We 
select all events that have two high- momentum tracks {\p\ > 70% -Ebcam) with- 
out associated high-energy showers (-Emax < 70% _Ebcam), that is, consistent with 
two high-energy muons accompanied by anything. We plot the invariant mass of 
these two muons in Figure 5.6, after subtracting continuum processes using the 



94 

off-resonance data. Muon pairs from direct T decays are easily distinguislied from 
cascade-to-leptons events. Tlie Monte Carlo reproduces the invariant mass spec- 
trum with only tiny errors in the calibration of the magnetic field (a horizontal 
shift in the plot). 

We measure Bcas relative to B^^ by fitting Monte Carlo cascade-to-muon pairs 
and Monte Carlo direct muon pairs to the data in Figure 5.6. This fit has two free 
parameters, the magnitude of i3cas and the magnitude of B^^. We assign conserva- 
tive 10% uncertainties to this procedure, which overwhelm i3^^ uncertainties and 
yield only 0.1% uncertainties in the final efficiency determination. The resulting 
-Bcas for T{2S) is (1.58 ± 0.16)% and for T{3S) is (1.34 ± 0.13)%, accounting for 
the factor of two from cascade-to-electron pairs. 

Applying these corrections to the T(2S') and T(3S') efficiencies, we obtain 
(96.18 toit ± 0.15)% and (96.41 to.tt ± 0.13)% for the T{2S) and T{3S), re- 
spectively. The first uncertainty is derived from the T(2S') -^ 7i^7i^T{lS) study 
and is common to all three resonances. The second uncertainty arises from our i3cas 
measurements and is independent for the T{2S) and T(3S'). The first uncertainty 
therefore cancels in ratios of Fee, while the second does not. 
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Figure 5.6: Invariant mass spectra of fi^fi~ in events with two energetic muons. 
Points with errorbars are continuum-subtracted data, and histograms are Monte 
Carlo simulations, with prompt T -^ fi'^fi^ shaded. The horizontal shift in energy 
scale is consistent with the uncertainty in our magnetic field estimate. 



Chapter 6 



Integrated Luminosity 

Now that we can reliably count the number of hadronic T decays at every Eqm we 
sampled, we need to determine the hadronic cross-section this implies by measuring 
the integrated luminosity of the same data. We get the integrated luminosity from 
a count of Bhabha events, since these are plentiful and their rate can be accurately 
calculated from perturbative QED. Integrated luminosity is the ratio of Bhabha 
events counted to the efficiency-weighted Bhabha cross-section, and the hadronic 
T cross-section is the ratio of hadronic T events to the integrated luminosity. 

Any theoretically calculable process can be used to determine the integrated lu- 
minosity; Bhabhas were chosen primarily for their abundance, since this minimizes 
statistical uncertainty. A Bhabha count is complicated by the fact that T -^ e+e~ 
is indistinguishable from Bhabhas on an event-by-event basis, and this background 
peaks under the resonance. Alternatively, one could determine the integrated lu- 
minosity from e+e^ -^ 77 because T y^ 77 (Electromagnetic decays do not violate 
parity). Unfortunately, this comes at the cost of a factor of 8.6 in statistics. The 
T -^ e+e^ contribution can easily be controlled, and the additional statistical 
power is valuable when determining ratios of ree{nS)/ree{iTLS), so we determine 
integrated luminosities from Bhabhas and use 77 events as a cross-check. 

6.1 Bhabha Count and 77 count 

To select Bhabha events, we require 

• two or more tracks with momenta between 50% and 110% of -Ebcam, and an 
energy sum (including showers from bremsstrahlung radiation as the elec- 
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trons propagate through the detector) of more than 90% Ecu- 

• The larger (smaller) track | cos^| must be less than 0.766 (0.707), and 

• each track must be associated with a calorimeter shower, with the larger 
(smaller) shower energy divided by track momentum (-Eshowcr/|Ptrack|) being 
greater than 80% (50%). 

With these cuts, backgrounds other than T -^ e^e~ are negligible. Different 
thresholds are set for the larger and smaller angles and calorimeter energies to 
reduce sensitivity to the threshold values, and possible efficiency variation with 
time. See Figure 6.f for an illustration of this kind of cut. 

Contamination from T -^ e^e^ at the resonance peaks is 3.8%, 1.4%, and 
1.0%, respectively. This background is readily calculated for any Eqm by multi- 
plying the T lineshape by Bee (we assume Bee = B^lfM for greater precision) and 
the cut efficiency for T -^ e^e^ (which has a different angular distribution than 
Bhabhas). Since the T lineshapes are derived from cross-section measurements, 
this is a circular dependence, so we applied an iterative procedure, starting with 
77 luminosity. The T and continuum e^e~ interfere, so we also calculate an in- 
terference term (Equation 2.8) with aint = 0.60, 0.87, and 0.69 for T(IS'), T(2S'), 
and T(3S'), respectively. The effective cross-section for e~^e~ as a function of Eqm, 
including the continuum and T contributions, is presented in Figure 6.2, with and 
without the interference term. The presence of the interference term has negligible 
impact on the lineshape fit results. 

To select 77 events for our cross-checks, we require 

• two showers (subscripted 1 and 2) with energies higher than 70% of -Ebeam, 

• I cot^i + cot ^2 1 < 0.1 (showers back-to-back in 6, the polar angle), and 



98 



04 




theta of shower/track 1 

Figure 6.1: The distribution of two anti-correlated variables with an asymmetric 
cut (the masked square regions). Horizontal and vertical smearing around the 
central diagonal is independent; we cut each track or shower in a way that doesn't 
imply a cut on the other. 
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Figure 6.2: The influence of T ^ e+e~ on effective e+e~ cross-section, with in- 
terference (solid) and without interference (dashed). Note that the vertical axis is 
zero-suppressed. 
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• I sin(0i — 02)1 < 0.04 (showers back-to-back in 0, the azimuthal angle). 

• The larger (smaller) | cot^| must be less than 1.28 (1.18) (this is within the 
calorimeter barrel), 

• the larger (smaller) | cot 6*1 must be greater than 0.15 (0.05) (avoiding the 
central region, where trigger efficiency is low), and 

• there must be no tracks in the event. 

When we select 77 events, we are dependent on only one trigger, barrel-bhabha, 
since this is the only trigger that doesn't require any tracks. We studied the effi- 
ciency of this trigger with Bhabhas, and found calorimeter tiles whose efficiencies 
dropped for a significant fraction of the data-taking period. Rather than apply- 
ing a run-dependent efficiency, we masked out these regions with our cuts. The 
largest shower on the western side of the detector must not be found in any of 
these regions: 

• -||7r < < ^TT and |cot^i + cot^2| < 1-08, 

• ~il^ < < -||7r and | cot^i + cot ^2! > 1-90, and 

• -0.4 < < -0.3. 

Given these angular cuts, the barrel-bhabha trigger is 99.67% efficient with only 
statistical deviations. (The efficiency of every run is above 99.2%.) 

6.2 Overall Luminosity Scale 

The efficiency-weighted Bhabha cross-section is the second ingredient in the lumi- 
nosity measurement. This sets the luminosity scale for all Bhabha counts. The 
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scale factor, a number of nb~^ per observed Bhabha, is the inverse of the efficiency- 
weighted Bhabha cross-section (nb). 

The efficiency-weighted cross-section is the cross-section of observed Bhabhas. 
Expressed as an integral over a single variable for clarity, 

a^i, = 1^ ^ e{9) d9 (6.1) 

where e(^) is our detector's Bhabha cut efficiency at a given polar angle 6. Rather 
than exhaustively simulating the detector's response to Bhabhas in 6 bins, we gen- 
erate Bhabhas with an angular cut-off beyond the detector's geometric acceptance 
(where e{9) = 0), calculate the cross-section this represents (ao), and multiply it by 
the efficiency of these simulated events (cq), determined by passing them through 
the detector simulation. Both o-q and Cq depend on our choice of cut-off, but 
the product doesn't. This product is the desired efficiency- weighted cross-section 
because 

.„ = / -de' and .„ = y -.(0)ie. (6.2) 

where dP/dO is the probability distribution of Bhabhas with polar angle 9 in our 
simulation, which is normalized: 

Equation 6.1 may be derived from Equations 6.2 and 6.3. 

We simulate Bhabhas with the Babayaga event generator, which calculates 
dP/dO and ctq to fourth order in the fine structure constant [17], with an angular 
cut-off of |cos^inax| = 0.819. Passiug these simulated events through our full 
detector simulation, we calculate an efficiency-weighted cross-section of 8.993 ± 
0.035 nb-i at an Ecm of 9.43 GeV, 7.945 ± 0.031 nb^^ at 10.00 GeV, and 7.361 
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± 0.032 nb-i at 10.33 GeV, which are the off-resonance T{1S), T{2S), and T{3S) 
energies, respectively. 

The Monte Carlo reproduces data distributions well at all three energies, which 
is demonstrated for the 10.33 GeV data in eight relevant distributions in Fig- 
ures 6.3-6.10. 

The efficiency-weighted Bhabha cross-sections at 9.43, 10.00, and 10.33 GeV 
differ by more than 1/s because the Monte Carlo finds these Bhabha cuts, partic- 
ularly the requirement that the energy sum of the two tracks be greater than 90% 
of Ecu, to be energy-dependent. We checked this claim by comparing Bhabha 
counts with 77 counts, and by loosening the energy sum cut. Both methods reveal 
the same 2% per GeV energy dependence from T(IS') to T(3S'). 

We additionally determine the efficiency-weighted cross-sections of e+e" -^ 
H^fi~ and e^e^ -^ 77 to reduce systematic uncertainties by comparing the lumi- 
nosity predicted for the same dataset by different processes. We follow the method 
of [18] to assign 1.6% systematic uncertainties for Bhabha efficiency-weighted cross- 
section, 1.6% for fi^fi^, and 1.8% for 77. The sources of these uncertainties are 
tabulated in Table 6.1: e~^e~ and fi~^fi~ uncertainties are dominated by the degree 
of resonance interference and the track-finding efficiency, while 77 uncertainties 
are dominated by the photon-finding efficiency and angular resolution. 

We can compare these as luminosity measurements by dividing the e~^e~ , /i^/i~, 
and 77 counts from the same sample — all available off-resonance data — by their 
efficiency-weighted cross-sections. We obtain the integrated luminosities plotted 
in Figure 6.11, which are in good agreement, considering their systematic uncer- 
tainties. The e^e~ and fi'^fi' modes share the same track-finding efficiency, so we 
draw error bars with this systematic removed in the Figure. 
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Figure 6.3: Largest track momentum divided by i?bcam in the 10.33 GeV data 
(points) and Monte Carlo (histogram) with other cuts apphed. Data between 0.5 
and 1.1 are accepted. 
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Figure 6.4: Second-largest track momentum divided by i?bcam in the 10.33 GeV 
data (points) and Monte Carlo (histogram) with other cuts applied. Data between 
0.5 and 1.1 are accepted. 
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Figure 6.5: Sum of two largest-momentum track energies and associated 
bremsstrahlung showers divided by Eqm in the 10.33 GeV data (points) and Monte 
Carlo (histogram) with other cuts applied. Data above 0.9 are accepted. 
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Figure 6.6: Positron cos^ (polar angle) distribution in the 10.33 GeV data (points) 
and Monte Carlo (histogram) with other cuts applied. The larger (smaller) | cos^| 
of the two electrons must be below 0.766 (0.707) for acceptance. 
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Figure 6.7: Largest -Eshowcr/|ptrack| divided by -Ebcam in the 10.33 GeV data (points) 
and Monte Carlo (histogram) with other cuts apphed. Data above 0.8 are accepted. 
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Figure 6.8: Second-largest -Eghowcr/IPtrackl divided by -Ebeam in the 10.33 GeV data 
(points) and Monte Carlo (histogram) with other cuts applied. Data above 0.5 are 
accepted. 
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Figure 6.9: Positron (azimuthal angle) distribution in the 10.33 GeV data 
(points) and Monte Carlo (histogram) with other cuts applied. This variable is 
not used for cuts. 
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Figure 6.10: Number of charged tracks in the 10.33 GeV data (points) and Monte 
Carlo (histogram) with other cuts applied. This variable is not used for cuts. 
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Table 6.1: Fractional systematic uncertainties in our determinations of the 
efficiency- weighted cross-sections of e+e~, f^~^fJ.~ , and 77. All values are percent- 
ages. 





e+e 


fi+fi- 


77 


Finite Monte Carlo sample 


0.4 


0.5 


0.6 


Radiative corrections 


0.5 


0.5 


0.5 


Resonance interference 


1.0 


1.0 




Trigger efficiency 


0.1 


0.1 


0.7 


Track-finding efficiency 


1.0 


1.0 




Photon-finding efficiency 






1.0 


Dependence on cuts 


0.5 


0.3 


1.0 


Cosmic ray backgrounds 




0.2 




ISR tail backgrounds 




0.1 




Total 


1.6 


1.6 


1.8 
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Figure 6.11: The integrated luminosity of all off-resonance data combined, as de- 
termined from e+e^, n^fi^, and 77 counts. Outermost errorbars include all sys- 
tematic uncertainties, and the second errorbars on e'^e~ and fi^fi" have common 
tracking and resonance interference systematics removed. Innermost errorbars (vis- 
ible only for fi^fj^^) are statistical-only. The weighted average and RMS of the three 
measurements are represented by dashed and dotted horizontal lines. 
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The weighted average luminosity from e^e^, /x^/x^, and 77 is only 0.2-0.6% 
higher than the e^e~ luminosity alone (depending on dataset), so we modify our 
luminosity determination to return the average luminosity from the three processes. 
The integrated luminosity is 

0.1114 (Ecm/9.43 GeV)^ nb^^ per observed Bhabha event, (6.4) 

0.1266 (i^civi/lO.OO GeV)^ nb~^ per observed Bhabha event, and (6.5) 
0.1361 (-Ecm/10.33 GeV)^ nb~^ per observed Bhabha event (6.6) 

in the three datasets. 

Without knowing that fi'^fi^ and 77 measurements reproduce the Bhabha re- 
sult, we would have a 1.6% uncertainty common to all three luminosity scale fac- 
tors. However, we can incorporate this information by assigning the e~^e~ , H~^H~ , 
and 77 RMS differences from the average as the common uncertainty. Thus, the 
luminosity scale factors have an uncertainty of 1.3%. The luminosity scale factor 
is a factor in reeFhad/rtot, so the fractional uncertainty in this scale factor adds to 
the Feerhad/rtot fractional uncertainty in quadrature. 

The systematic uncertainties in the three luminosity scale factors are not inde- 
pendently 1.3%: they will partly cancel in ratios. We can see this by considering 
the minimum information necessary to determine rge{ns)/rge{iTiS). One must 
know the ratio of Bhabha cut efficiencies near the T[nS) and the T(m5') and the 
scaling of Bhabha cut efficiency with Ecu, but the conversion from the number of 
Bhabha events to inverse nanobarns will cancel. Thus, the dominant uncertainty in 
ratios of luminosity scale factors is the ratio of Bhabha cut efficiencies, determined 
with 0.5% uncertainty by Monte Carlo. 

To blind our analysis, we determined the overall luminoisty scale factors last. 
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after all background (Chapter 4), efficiency (Chapter 5), and beam energy studies 
(Chapter 7) were completed. Until that point, we only knew the ratios of luminosi- 
ties of our data samples, and therefore only T ee{nS) /T eeimS) . We incorporated 
these scale factors (using 77 event counts, rather than Bhabha event counts) just 
before presenting preliminary Fee results at the European Physical Society meeting 
in July, 2005. 

6.3 Consistency of Bhabhas with 77 

Above, we took advantage of the 77 luminosity measurement's very different sys- 
tematic uncertainties to check and correct the Bhabha scale factors off-resonance, 
but we have not yet used the fact that 77 counts are unaffected by T decays. In 
this Section, we will compare Bhabha and 77 rates as a function of Eqm through 
the T resonances, to test the T -^ e+e~ correction. We tune the 77 luminos- 
ity measurement to yield the same results as Bhabhas off-resonance and measure 
the 77 and Bhabha luminosities in T data. Both methods should yield the same 
integrated luminosity, within statistical uncertainties. 

In Figure 6.12, we plot ratios of luminosity measurements for off-resonance, 
scan, and peak data on the T(IS'), T{2S), and T(3S'), and observe a discrepancy 
on the T(IS') and T(3S'). The 77 luminosity scale factor is tuned to reproduce 
Bhabha luminosity off-resonance, but not separately for the three off-resonance 
datasets, so the discrepancy we observe occurs only at the T resonances. It is not, 
therefore, linear with respect to Eqm, as would be expected if the discrepancy were 
due to energy-dependent cut efficiencies. (The discrepancy is also 20 times larger 
than the energy dependence in the Bhabha cuts, and the 77 cuts have no energy 
dependence on this scale.) 
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Figure 6.12: Ratio of Bhabha luminosity to 77 luminosity, with the weighted 
average continuum ratio set to unity. 



Ill 

The direction of this effect, and the fact that it only influences T data, could 
be explained if we were over-subtracting T -^ e+e^ from our Bhabha count. How- 
ever, the relative magnitudes of the T(IS') and the T(3S') discrepancies cannot be 
accounted for. An over-subtraction error would be proportional to peak hadronic 
cross-section times Bee, which is 0.45 nb for the T(IS') and 0.096 nb for the T(35'), 
but the discrepancies appear to be equal. Even considering the uncertainties in 
these measurements, the over-subtraction hypothesis is ruled out by 8.7 standard 
deviations. 

We do not know whether the discrepancy is due to an error in the Bhabha 
measurement or in the 77 measurement, so we apply half of this discrepancy as 
a correction and add half the discrepancy and its uncertainty in quadrature to 
the total uncertainty in luminosity, to cover the ambiguity. This is only a 0.4% 
uncertainty in the luminosity of each resonance, and therefore a small contributor 
to the total Fee uncertainty. 



Chapter 7 

Beam Energy Measurement 

Having presented all the details necessary to properly measure the hadronic T 
cross-section, the vertical axis in our lineshape plot (Figure 2.2), we now consider 
our measurement of Eqm, the horizontal axis. In Chapter 3, we discussed the 
mechanism which determines the beam energy of each run. While this reckoning 
differs from the true beam energy by 18 MeV near the T resonances, we will show 
in this Chapter that measurements of beam energy differences are robust enough 
for the precision demands of this analysis. 

The masses of the T resonances have been measured with 0.3-0.5 MeV precision 
at Novosibirsk [19], so we can use the T lineshapes as calibrating markers in beam 
energy. If we correct our Eqm measurements by a constant shift at the T(IS'), we 
almost reproduce the Novosibirsk T(2S') and T(3S') masses: our average T(2S') 
mass is 0.65 ± 0.18 MeV too low, and our average T(3S') mass is 1.60 ± 0.30 MeV 
too low (Figure 7.1). If this error is purely a function of beam energy and is 
not related to alterations in CESR's configuration between the T(3S'), T(IS'), 
and T(2S') data-taking periods, then our beam energy measurements differ from 
the true beam energies by a linear transformation, which we determine with a 
fit to the three apparent T masses in Figure 7.1. If this tranformation applies 
to small changes in beam energy (scan point differences, about 10 MeV) as well 
as large (300 MeV between resonances), the resonance scans are too narrow, and 
Teerhad/rtot IS too low, by ouly 0.2%. Because of the uncertainties noted above, 
we do not apply this as a correction. 

Potentially more worrisome are variations in the beam energy with time. The 
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Figure 7.1: Difference between the measured mass and the true mass of the T(IS'), 
T{2S), and T(3S'), where the T(IS') measurement has been shifted to the true 
mass. Outer errorbars represent the RMS of all scan measurements at each reso- 
nance, and inner errorbars are statistical-only. The dashed curve is a linear fit. 



114 

measurement of beam energy differences is sensitive to the placement of the NMR 
probes in the test magnets, because at the 0.05% level (corresponding to 1 MeV in 
Ecu), the magnetic field in the test magnets is a function of position. Every week, 
these probes are exposed to the possibility of movement as a consequence of CESR 
machine studies, since these studies may involve reconfiguring the test magnet. 
To protect our lineshape scans from discrete shifts in beam energy calibration, we 
divided our data-taking into short, independent 10- hour scans, plus 38 hours of 
subsequent peak running, separated by about a week. The test magnets were not 
disturbed during these dedicated scans or during the peak data-taking associated 
with each scan. Changes in the apparent T mass from one week to the next alert 
us to shifts in the beam energy measurement between scans, which we observe as 
a slow drift on the order of 0.5 MeV per month (Figure 7.2). To be insensitive 
to these drifts, we include the apparent T mass of each 10+38-hour scan as an 
independent parameter in the lineshape fit. Values of Ecu in the lineshape fit are 
only relative to the apparent T mass from the week the data were taken. 

We are only sensitive to random changes in the beam energy calibration (true 
-E-beam — mcasured -Ebcam), or jitter, on a 10-hour timescale. We have already en- 
sured our insensitivity to week-by-week fluctuations. If the beam energy fluctuates 
on very short timescales, much less than an hour (the length of a run), then the 
jitter only contributes to beam energy spread. (In this sense, the beam energy 
fluctuates by about 4 MeV with every collision!) We reduce our sensitivity to 
monotonic 10-hour calibration drifts by alternating scan points above and below 
the T peak, so that a drift would not systematically widen or narrow the lineshape. 

To check for changes in the beam energy calibration on the 10-hour timescale, 
we measure the cross-section at a high- derivative point on the lineshape twice, usu- 
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Figure 7.2: Difference between the measured mass and the true mass of the T{1S), 
T{2S), and the T(3S') (top to bottom), for each individual scan. Dashed horizontal 
lines represent the weighted averages. 
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ally at the beginning and end of a scan. Since the derivative is high, significantly 
different cross-section results at the same reported beam energy would be an in- 
dication that the actual beam energy has shifted between the two measurements. 
Any cross-section difference (o"i — (T2) can be converted into a true beam energy 
difference {Ei — E2) at a given lineshape derivative [da/dE) by 

To find possible shifts in calibration, we subtract the true beam energy difference 
from the reported difference, and propagate uncertainties from the cross-section 
measurement. We call this a shift (sj ± 6sJ- The shifts (plotted in Figure 7.3) are 
not quite consistent with zero, as the total x^ = Xlil-^i/^sJ^ is 31.2 for 29 degrees 
of freedom (a 0.64% confidence level). This broadening of the {sj} distribution 
can be accomodated by a uniformly-random jitter j. To quantify j, we construct a 
log-likelihood function of the {sj} data, assuming them to be Gaussian-distributed 
with statistical and jitter components to the width. 

Only jitters larger than 0.05 MeV are inconsistent with the data at the 68% con- 
fidence level (the value of j necessary to lower the log-likelihood by 1/2; see Fig- 
ure 7.4). Alternatively, we could have constructed an S*- factor in analogy with the 
Particle Data Group's method of calculating uncertainty from a set of mutually- 
inconsistent experiments: 

In this formulation, the maximally-allowed jitter is the value which reduces S{j) 
to unity, which is j = 0.07 MeV. Both methods yield roughly the same value. 
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so we can be confident that the beam energy cahbration jitter is approximately 
0.07 MeV, which is 7 parts per milhon. 

To learn what fluctuations of this size imply for uncertainty in Fgerhad/rtot, 
we simulated the lineshape fits with random perturbations in the measured beam 
energy. The nominal -Ecm and luminosity of each simulated cross-section mea- 
surement were copied from the real data sample, but the cross-sections themselves 
were derived from an ideal curve with statistical errors. Without perturbing the 
simulated beam energy measurements, the fits return our input Feerhad/rtot with 
perfectly-distributed fit x^ values. When Eqm is randomly perturbed with a stan- 
dard deviation of 0.07 MeV, the value of TeeFhad/rtot fiuctuates up and down by 
0.2% of itself, and the x^ increases by only 5-30 units. We therefore assign a 0.2% 
systematic uncertainty in reeFhad/rtot and Fee due to beam energy measurement 
errors. 
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Figure 7.3: Beam energy shifts (sj) determined by pairs of repeated cross-section 
measurements, plotted with respect to the time between the first and the second 
measurement. The weighted mean of {si} is 0.02 ± 0.03 MeV, so the apparent 
vertical asymmetry of this distribution is an illusion. 
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Chapter 8 



Lineshape Fitting 

We have reached the core of the analysis, the fits of the T hneshapes. All previous 
studies either deliver data to these fits or define the curve that the data are fit to. 
In this Chapter, we will review the fit function and all of its parameters, present 
the fit results, and discuss hadron-level interference as a fitting issue. 

8.1 The Fit Function 

8.1.1 Hadronic Peak 

The central feature of the fit function is a Breit-Wigner curve representing the 
hadronic cross-section (Equation 2.3), convoluted with beam energy spread and 
initial-state radiation (ISR). The beam energy spread is modeled by a Gaussian 
with unit area, and the ISR distribution is calculated to fourth order in perturba- 
tive QED by Kuraev and Fadin [20] (Equation (28)). 

Most of the floating parameters of the fit modify this contribution. The area 
of the Breit-Wigner is the most important floating parameter, as this is how we 
determine reeEhad/rtot and Fee (Equation 2.2). To let each 48-hour scan's T mass 
float independently (Section 7 on page 114), we fix the mass in the fit function 
and allow the Eqm measurements in the data to be shifted scan-by-scan. This 
way, when we plot the lineshape fits, a single curve represents the fit to all data. 
In Chapter 3, we discussed the variability of beam energy spread (Section 3.1 on 
page 31). To allow scans to have different beam energy spreads (Gaussian widths), 
we use different fit curves. We do not transform the data in analogy with the mass 
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shifts because the ISR tail spoils the linearity of the transformation. 

The full width F is a parameter in the Breit-Wigner distribution, but we do 
not allow this parameter to float. Our fitted Feerhad/rtot is sensitive to F only at 
the 0.03% level, so we fix each F to its previously measured value [21]. 

The hadronic T efficiency is another multiplicative constant in our fit function. 
We multiply our fit function by efficiency rather than dividing our data by effi- 
ciency because our data includes several different components, each with a different 
efficiency. 

We add interference with the continuum (Section 2.1 on page 22) to the Breit- 
Wigner before convolution. That is, the signal lineshape is 

^res+i„t(^CM) = Ures(^') + &int{E')] ^G^ISR (8.1) 



where cTi-cs and 5"int are defined in Equation 2.5 and represents convolution with 
the beam energy spread Gaussian (G) and the ISR distribution (ISR). This (Xrcs+int 
does not represent a physical cross-section until we add the continuum piece. (It 
is negative for some values of Ecu-) The interference between resonant and con- 
tinuum qq decays is characterized by two constants, 0o and aint (Equation 2.8), 
both of which are known and do not need to fioat in the fit. 

8.1.2 Tau-Pair Peak 

The T~^T^ background term has the same form as the hadronic signal term: it is 
the sum of a Breit-Wigner and 5";^^^ , convoluted by the same beam energy spread 
and ISR distribution. This term introduces no new fioating parameters: the Breit- 
Wigner area for t'^t' is the Breit-Wigner area for hadronic T, multiplied by a 
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ratio of branching fractions and efficiencies 

= -j^ X («.2j 

areahad iJh&d Chad . 

The value of Err we use in the fit function is from [22]. We use this value to 
determine reeFhad/rtot, but when we determine Fee, we perform a separate fit with 
Brr = B^^ and B^^ from [12]. This allows us to subtract t^t^ backgrounds with 
the same branching fraction that we use to multiply it back in, when we convert 
reeFhad/rtot to Fge (Equatiou 2.10). 

The magnitude of interference, aJnt^ 5 is also different for the r^r^ peak, be- 
cause a{e^e^ -^ r^r^) and F(T -^ t^t^) are not equal to (j{e^e^ -^ qq) and 
F(T -^ qq), respectively. 

8.1.3 Background Terms 

We add a single 1/s term to the fit function to represent continuum qq, radiative 
Bhabhas, r^r^ , and any residual Bhabha or /x^/i" backgrounds. The magnitude 
of this term fioats in the fit, though it is primarily infiuenced by our large off- 
resonance data point. 

The resonance interferes with only part of the contniuum, so the magnitude of 
resonance interference is not tied to the fitted value of the 1/s term. The qq and 
T^r^ cross-sections that enter into ctint and aj^t^ are either derived from other 
experiments ([9] for a[^t) or are calculated from QED (for aj^t^ )• 

This dissociation of the continuum cross-section from the interference term is 
artificial, since interference is only meaningful in the presence of both resonance 
and continuum amplitudes. It would be more natural, for instance, to include the 
continuum cross-section in the beam energy spread and ISR convolutions with the 
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resonance and interference terms, since the natural lineshape includes all three 
terms before it is smeared by the physical beams. However, the continuum cross- 
section doesn't depend on Eqm sharply enough for this to matter. 

For T(2S') and T{3S) fits, we add the fitted lineshape(s) from T(IS') and T(2S') 
as backgrounds, because the ISR tails from these resonances introduce 0.4-0.6% 
corrections. These background components do not fioat in the fit. (We fit the 
three resonances separately.) 

The log s correction motivated by two-photon fusion backgrounds is tied to the 
overall continuum normalization through / = 0.080 ± 0.005, the fraction of logs 
component to 1/s at 9 GeV, determined by a fit from Section 4.3.1 on page 68. 
The total continuum cross-section is therefore 

.»(^CM) ^ .».. (d - /)^ . f log ^) (8.3) 

where 0"°^^^ is the fioating parameter. 

We allow cr^ont ^o float independently in the T(IS'), T(2S'), and T(3S') fits. 
Since it represents the cross-section at a fixed Eqm, the three fits ought to return 
the same value. Instead, the T(IS'), T(2S'), and T(3S') fits yield 9.355 ± 0.010 nb, 
9.318 ± 0.007 nb, and 9.315 ± 0.011 nb, respectively. The T{2S) and T{3S) values 
are consistent with each other, but the T(IS') fit has a 0"°^^^ which is 3.0 standard 
deviations higher than T(2S'). However, this difference in cr'^^^^ values is only 0.04% 
of o"J?o^^ itself. It is possible that our ISR tail corrections are too large, an error 
which is absorbed into o"°o„t, rather than biasing the fit for Tgerhad/rtot- 

8.1.4 Summary of all Floating Parameters 

The T(IS') fit has a total of sixteen parameters: reeFhad/rtot, Mi, . . . Mu for its 
eleven 48-hour scans, AEi, AE2, and AE3 for its three groups of scans with poten- 
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tially different beam energy spreads (see Table 3.1), and aj^o^^ for the continuum 
background normalization at 9 GeV. 

The T(2S') fit has a total of nine fioating parameters: Feerhad/rtot, ^i, • • • , Mq 
for its six 48-hour scans, a single beam energy spread AE (the CESR orbits were 



stable), and a^ 



cont ' 



The T(3S') fit has a total of sixteen fioating parameters: Teerhad/rtot) Mi, . . . , Mj 
for its seven 48-hour scans, AEi, . . . ,AEi, a beam energy spread for each scan. 



and a°„nt. 



8.2 Systematic Uncertainties from Constants in the Fit 

Each constant in the fit function may introduce systematic error, but in every 
case, this error is negligible. Kuraev and Fadin estimate that their ISR calculation 
has a 0.1% uncertainty; varying the tail normalization by this amount changes our 
FeeFhad/rtot rcsult by 0.05%. As previously mentioned, the uncertainty in F affects 
our result by less than 0.03%. Uncertainty in the qq interference magnitude (from 
R) yields a 0.02% uncertainty, and the t~^t~ interference magnitude (from Brr) 
yields 0.08%. Uncertainty in the two-photon fraction / yields 0.002%. The sum of 
all of these uncertainties in quadrature is 0.10%. 

We do not propagate the uncertainty in the t~^t~ efficiency (first mentioned in 
Section 8.2 on page 123) because it is multiplied by Brr or B^^ (in the FgeFhad/Ftot 
and Fee fits, respectively). The fractional uncertainty in the Monte Carlo modeling 
of T^T^ efficiency is not significantly greater than 2%, the smallest Brr or i3^^ 
fractional uncertainty. 

The hadronic efficiency also enters the fit function as a constant, but we will 
treat this as a major systematic uncertainty in the next Chapter. 
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Table 8.1: Best-fit Breit-Wigner area and x^ significance. 



J a{e^e ^ T ^ hadronic) dii^cM reduced x^ confidence level 



T{1S) 324.12 ± 0.92 MeV nb 

T{2S) 133.65 ± 0.94 MeV nb 

T{3S) 89.25 ± 0.88 MeV nb 



240.4 



203 - 16 


= 1.29 


0.51% 


107.2 
75-9 


= 1.62 


0.10% 


154.5 
175 - 16 


= 0.97 


58.6% 



8.3 Fit Results 

We used a C++ version (1.5.2) of MINUIT from the SEAL package [23] to perform 
a x^ fit of our cross-section data to our fit function. The resulting best-fit values 
for the Breit-Wigner area used in determining Fee are presented in Table 8.1. The 
best-fit functions are plotted in Figure 8.1, and Figures 8.2-8.4, overlaying data. 
At the top of Figure 8.1 and in Figures 8.5, 8.6, and 8.7 are pull distributions. 
Pull is the ratio of residual to uncertainty, expressing the significance of deviations 
from the best-fit line. All data points are visible in these plots; no point has a 
pull greater than 3.5. Best-fit T masses and beam energy spreads are plotted in 
Figures 7.2 and 3.2, respectively. 

This fit function describes our data well: no trends are evident in pull versus 
energy, but the T(IS') and especially the T(2S') pull distributions are wider than 
would be expected from statistical fiuctuations. We can also see this by noting 
that the x^ per degree of freedom (A^dof) is improbably high (Table 8.1). In the 
T(2S') data, the high-energy tail is the largest deviation, contributing 9.1 units to 
the total x^) but this is only a fifth of the excess. (Dropping this point changes 
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Figure 8.1: Best-fit lineshapes (solid) to the T(IS'), T(2S), and T{3S) data (points), from which Teerhad/rtot is determined. 
Dashed curves are total background estimates, and the pull distributions above each plot indicate the statistical significance 
of all deviations from the best-fit line. Insets for each resonance present a tail measurement 100 MeV, 60 MeV, and 45 MeV 
above the T(IS'), T(2S'), and T(35') masses. (Fit functions with multiple beam energy spreads are represented by their 
average beam energy spread.) 
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Figure 8.2: Best-fit T(IS') lineshape, overlaid upon data from individual scans. 
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Figure 8.3: Best-fit T(2S') lineshape, overlaid upon data from individual scans. 
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Figure 8.4: Best-fit T(3S') lineshape, overlaid upon data from individual scans. 
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Figure 8.5: The pull distribution of T(IS') linesliape fits (a) as a function of Ecu, 
(b) as a function of date, and (c), as a histogram, fitted to a Gaussian curve, which 
is almost two standard deviations wider than unity. 
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Figure 8.6: The pull distribution of T{2S) linesliape fits (a) as a function of i?cM, 
(b) as a function of date, and (c), as a histogram, fitted to a Gaussian curve, which 
is almost two standard deviations wider than unity. 
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Figure 8.7: The pull distribution of T{3S) linesliape fits (a) as a function of Ecu, 
(b) as a function of date, and (c), as a histogram, fitted to a Gaussian curve, which 
is statistically consistent with unity. 
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the T(2S') Feerhad/rtot by 0.4%.) The simulations of fits with jittering beam 
energy, described at the end of Section 7 on page 117 demonstrate that uncertain 
beam energy measurements can increase the fit x^ by 5-30 units. If the T(IS') 
and T(2S') x^ values have both been artificially raised 30 units by jittering beam 
energy measurements, their natural confidence levels would be 10-15%, but there 
is no conclusive evidence that beam energy jitter is the cause, and assuming an 
increase of 30 units of x^ is the most extreme hypothesis. 

The larger-than-expected deviations are symmetric around the best-fit line for 
all Ecu and dates, and they do not favor a particular Ecu or date (Figures 8.5- 
8.7), so we treat them as though we had underestimated our statistical uncertainty 



(o"stat)- We add o-gtatv X^/^dof — 1 to our systematic uncertainty in quadrature, 
which has the same effect on the total uncertainty as if we had multiplied o"stat by 

Note that our fit does not suggest that the beam energy spread distribution 
is distorted from a pure Gaussian. If this were the case, we would observe trends 
in pull versus Ecu- Therefore, our assumption of a Gaussian beam energy spread 
is valid for our level of precision (statistical plus the above-mentioned systematic 
uncertainty, which is roughly ±0.1 nb). 

8.4 Hadron-Level Interference 

In our fit function, we assumed that continuum and T decays only interfere at the 
parton level, and not at the hadron level. That is, we assume that e+e^ -^ qq can 
interfere with e+e^ ^ T ^ qq but not e+e^ ^ T ^ ggg, even though both qq 
and ggg hadronize into some of the same final states. 

This gg-only scheme optimizes the T(IS') fit, assuming no phase difference 
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between resonance and continuum at Ecu ^ Mr (0o = 0). We determined this 
by repeatedly fitting our lineshape, assuming different interference magnitudes 
ctint (defined in Equation 2.8). We plot the x^ oi these fits in Figure 8.8. The 
interference magnitude ctint is proportional to the square root of the branching 



fraction of T to the interfering final state, which is y/Bi^t = V8.9% for qq. Our fit 
is optimized by ajnt = 0.016 ± 0.004, or Bi^t = 7.8 ± 4.2. The x^ of the T{1S) fit 
with no interference is 15 units larger than our gg-only scheme, which corresponds 
to 3.7 standard deviations. 

However, if e^e^ ^ T — ;► ggg — > hadrons and e~^e~ —* qq —^ hadrons signifi- 
cantly interfere, the phase difference 0o need not be zero. Strong and Electromag- 
netic T decays to hadrons can differ in phase, an effect which has been observed 
in exclusive decays of charmonium [24] [25]. This Strong minus Electromagnetic 
{ggg — qq) phase difference is constant with respect to Ecu, which gives the res- 
onance minus continuum phase difference a constant offset as it evolves through 
the resonance peak. In this case, 0o need not be zero, and our fit is less well 
constrained. 

The question of whether ggg and qq interfere inclusively (in the sum over all 
hadronic final states) is fundamental and interesting in its own right, concerning the 
quantum mechanics of the hadronization process. Unfortunately, our data cannot 
answer this question for all phases, because if 0o ~ i^/2, the interference term a^it 
has the same Ecu dependence as the resonance Breit-Wigner (Equation 2.8). It has 
been argued [26] that inclusive ggg/qq interference is a small effect because every 
exclusive final state may have a different ggg/qq phase. The sum of interference 
cross-terms includes cancellations and therefore grows more slowly with the number 
of exclusive final states than the direct cross-section terms. 
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Figure 8.8: A fit of the T(IS') for the magnitude of hadronic interference (aint), 
assuming 0o = 0. The points are T(IS') fits with trial ajnt values, the vertical axis 
is the x^ of these fits, and the solid line is a parabola drawn through the points. 
Dashed and dotted lines indicate the minimum and 68% confidence level bands, 
and the circled point represents the gg-only fit we used in this anaysis. 
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Though improbable, it is possible that exclusive final states add coherently, such 
that the inclusive ggg and qq cross-sections significantly interfere with ggg — qq 
phase angle (pggg- The magnitude of this interference cannot be maximal, since 
some ggg final states are distinguishable from qq decays. We define fggg as the 
fraction of the ggg amplitude which interferes with qq. If (f)ggg = 0, our previous 
study applies (Figure 8.8) and fggg < 0.12 at 95% confidence level. We can also 
put limits on fggg which are independent of (pggg by considering the fact that ggg 
and qq have different quantum numbers. One-third of the qq amplitude is cc, but 
T(IS') — > ggg -^ cc is negligible [27]. Moreover, Strong decays preserve the T 
meson's isospin of / = 0, while Electromagnetic decays through a virtual photon 
have / = or 1. At most half of the uu and dd has /^ = 0, and therefore overlap 
with ggg. Adding these fractions yields 

/ 



cc ss 



1x0 + i xl 



dd uu 



1 I 2 



|xu+ |Xi + \| + |;x0.5 = 42% 

or fggg ^ 0-4- 

Our T(IS') fits constrain fggg as a function of (pggg. Assuming (pggg values 
between and 2tt, we fit for the maximum fggg consistent with the data at the 
68% confidence level (the value which raises the fit x^ by one unit), and plot 
this allowed region in Figure 8.9. For most of the (pggg range, our fit is more 
restrictive than the constraint we derived from fiavor and isospin considerations, 
though within 0.2 rad of ±7r/2, our fit yields no information. Outside of this range, 
the minimum x^ is attained within < fggg < 0.1. 

Hadron-level interference also affects our FgeFhad/rtot result, especially for (pggg 
near ±7r/2. From the same fits, we derive the following relationship between the 



136 
correction that must be applied to Feerhad/rtot and fggg and (pggg. 

A(reerhad/rt0t) = -5.42% T/^,, Sm{(Pggg + 0.3y\ +0.U%(fggg + 0.92 fgggA (8.5) 

Applying the limits on fggg shown in Figure 8.9 to Equation 8.5, we obtain the 
68% confidence level upper limits on T(IS') FeeFhad/rtot corrections shown in Fig- 
ure 8.10. Corrections for the T{2S) are 70% as large as this and corrections for 
the T{3S) are 65% as large. 

Complete a priori ignorance of inclusive ggg/qq interference would imply a 
2.0% uncertainty in FeeFhad/rtot and Fee for the T(IS'), a 1.4% uncertainty for the 
T(2S'), and a 1.3% uncertainty for the T(3S'). However, the possibility that fggg is 
as large as 0.4 is not widely believed to be likely, nor is it clear that (pggg ^ ±7r/2 
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Figure 8.9: The fraction of ggg amplitude allowed to interfere with continuum qq 
ifggg), accordlug to the T(IS') fit at 68% confidence level. We have no sensitivity 
to this parameter if the ggg — qq phase difference in T decays {4>ggg) is near ±7r/2. 




Figure 8.10: The possible correction to FeeFhad/rtot implied by our ignorance of 
fggg d^t 68% coufideuce level (gray) as a function oi(f)ggg. The dashed lines represent 
upper limits including isospin and fiavor arguments. 



Chapter 9 



Results and Conclusions 

In this Chapter, we present final results for Teerhad/rtot, for Fee, for T ee{nS) /T ee{rnS) , 
for r, and for \ip{0, 0, 0)p, all of which are derived from Feerhad/rtot- But first, we 
review the systematic uncertainties in FgeFhad/rtot- 

The hadronic T efficiency and integrated luminosity scale factor together deter- 
mine the multiplicative scale for cross-section, so fractional uncertainties in these 
two factors add in quadrature to the fractional uncertainty in FeeFhad/rtot- Un- 
certainties in backgrounds contribute negligibly to systematic uncertainty (their 
largest effect is statistical). At the end of Chapter 7, we determined the effect 
of beam energy uncertainty on the lineshape fits through simulations, and in Sec- 
tion 8.2, we determined the uncertainty due to the parameterization of our fit 
function. Our data have two unexplained features: Bhabhas and 77 events do 
not predict the same luminosity as a function of Ecu, and the fit x^ is signifi- 
cantly higher than the number of degrees of freedom for the T(IS') and T{2S). 
We have quantified both of these as systematic uncertainties. All uncertainties in 
Teerhad/rtot and Fee are listed in Table 9.1, including uncertainty from the correc- 
tion for leptonic modes (Equation 2.10), which applies to Fee and not FeeFhad/rtot- 

Our values of FeeFhad/rtot? quoted with statistical uncertainties first and sys- 
tematic uncertainties second, are 

(FeeFhad/rtot)(15) = 1.252 ± 0.004 ± 0.019 keV, (9.1) 

(reerhad/rtot)(25) = 0.581 ± 0.004 ± 0.009 keV, and (9.2) 

(FeeFhad/rtot) (35) = 0.413 ± 0.004 ± 0.006 kcV. (9.3) 
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Table 9.1: All uncertainties in rgeFhad/rtot and Fee- The "correction for leptonic 
modes" applies to Fee only, and "common hadronic efficiency" and "overall lumi- 
nosity scale" are common to all three resonances. 

Contribution to Fee T(15) T{2S) T{3S) 

Common hadronic efficiency (Section 5.2) 

T(2S'), T(3S') efficiency corrections (Section 5.3) 

Overall luminosity scale (Section 6.2) 

Bhabha/77 inconsistency (Section 6.3) 

Beam energy measurement drift (Chapter 7) 

Fit function shape (Subsection 4.3.1 and Section 8.1) 

X^ inconsistency (Section 8.3) 

Correction for leptonic modes (Equation 2.10) 

Total systematic uncertainty 1.5% 1.6% 1.5% 

Statistical uncertainty (Table 8.1) 0.3% 0.7% 1.0% 

Total 1.5% 1.8% 1.8% 
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Correcting for leptonic modes with a factor of (1 — 3i3^^), we obtain Fee: 

iTee)ilS) = 1.354 ± 0.004 ± 0.020 keV, (9.4) 

(ree)(25) = 0.619 ± 0.004 ± 0.010 keV, and (9.5) 

(Fee) (35) = 0.446 ± 0.004 ± 0.007 keV. (9.6) 

The total uncertainties for each resonance is less than 2%. Figure 9.1 compares 
our FeeFhad/rtot to all previous measurements: we find it to be consistent with, 
but more precise than, the world average. 

As a consequence of the measurement technique, our Fgg values represent the 
rate of T -^ e+e^ decays with no photons in the final state, not even very soft 
photons. This is because final-state radiation in Fee (TT — > 7e+e~) corresponds to 
initial-state radiation in the measured cross-section (e+e^ -^ l'^), which we ex- 
cluded with our fitting technique. Our result does include Electromagnetic vacuum 
polarization, a few-percent effect in which the virtual photon connecting e"'"e~ to 
T is interrupted by fermion loops, because this affects the decay process and the 
production process equally, and we never made any correction to remove it. 

This analysis represents a substantial gain in precision, but it is also more 
general than previous analyses in that we do not assume an T decay model to 
determine our efficiency. In previous analyses, efficiency was determined by se- 
lecting events from a Monte Carlo simulation, so potential errors due to modeling 
hadronization, unknown decay modes, or the detector simulation were not rep- 
resented by their quoted efficiency uncertainties. We relax this assumption for 
T(IS') decays and only assume that our Monte Carlo correctly scales from T(IS') 
to T(2S') and T(3S'). In addition, we obtain an upper limit on all-neutral T(IS') 
decays. The T(IS') branching fraction to events that generate zero tracks with 
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Figure 9.1: Comparison of our Feerhad/rtot with previous measurements. The dashed and dotted hues are world averages 
and uncertainties, excluding our new measurements [21]. 
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I cos^l < 0.93 and p± > 150 MeV is less than 1.01% at 90% confidence level. See 
Section 5.2.1 on page 85 for more details. 

This is the first analysis with sufficient precision to observe interference between 
the T(IS') resonance and the continuum. This effect was observed for the J/ip 
and ip' mesons in the fi'^fi~ channel within a year of their discovery [28] [29] 
because these charmonium resonances have much larger cross-sections relative to 
continuum. We therefore had the first opportunity to explore the nature of this 
interference in the T meson. We found evidence of and corrected for interference 
in the qq final state (Figure 8.8) and observed interference in the n^fi^ final state 
(Figure 2.4), but not in inclusive hadrons. 

These measurements provide useful constraints on the bb potential or checks 
on the validity of the approximation [30]. The degrees of freedom in these cal- 
culations are the T masses and widths. Our Fee measurements compliment the 
Novosibirsk measurements of the T masses [19], forming a set of six experimental 
inputs with better than 2% precision each. The T(3S') measurement is partic- 
ularly valuable, as it is the second ever performed (Figure 9.1), increasing the 
world knowledge by a factor of five. Also, the T(3S') wavefunction probes more of 
the non-Coulomb "confinement" potential than the T(IS') or T(2S'). The T(4S') 
wavefunction reaches further, but it is more complicated to extract Fge due to its 
vicinity to the B^B^ threshold [31]. The full width of the T(4S') is a function of 
Eqm and is no n- negligible (F ^ 20 at the T(4S') mass), making it necessary to 
model a QCD process to express this function. 

Finally, the T di-electron widths are valuable as a test of Lattice QCD. The 
ratios T ee{nS) /T ee{fnS) are especially useful, since the renormalization factor for 
the virtual photon current, -^match cancels, which allows us to test this aspect of 
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the calculation in isolation. Conversely, the ratios test the rest of the calculation, 
particularly the NRQCD treatment of b quarks and the staggered- quark formalism 
for virtual light quarks, in isolation and with higher precision than what can be 
attained with the absolute Fee calculations. These aspects are shared with the 
calculation of j'b, which is used to extract Vtd from B^-B^ mixing. 
Our experimental values of ree(nS')/ree(^5') are 

ree(2S)/ree(15) = 0.457 ± 0.004 ± 0.004, (9.7) 

ree(35)/ree(15) = 0.329 ± 0.003 ± 0.003, and (9.8) 

Tee{'iS)/Vee{2S) = 0.720 ± 0.009 ± 0.007, (9.9) 

where we have canceled systematic uncertainties which are shared among the pairs 
of resonances. These shared systematics are the common hadronic efficiency and 
the overall luminosity scale in Table 9.1. We overlay our ree(2S')/ree(15') experi- 
mental result on the HPQCD calculation (Equation 1.2) in Figure 9.2. 

We now use our Fee measurements to determine basic parameters of the T 
mesons: their total decay rates (lifetimes) and wavefunctions at the origin (sizes). 
The only known experimental access to these parameters is through Fee- Again 
assuming Bee = Bfj,f^, we find 

T{IS) = 54.4 ± 0.2 ± 0.8 ± 1.6 keV, (9.10) 

T{2S) = 30.5 ± 0.2 ± 0.5 ± 1.3 keV, and (9.11) 

T{3S) = 18.6 ± 0.2 ± 0.3 ± 0.9 keV. (9.12) 

The first two uncertainties are statistical and systematic, and the third uncertainty 
is propagated from B^^. The T(2S') and T(3S') values are lower than the averages 
quoted in [21] because the 2005 measurements of B^^ [12] are higher than previous 
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Figure 9.2: HPQCD calculations of ree(25')/ree(lS') times the ratio of masses 
squared as a function of lattice grid size squared. The extrapolation to zero lattice 
size and our new measurement are overlaid. 
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measurements and were not included in [21]. The above widths correspond to the 
following lifetimes: 

t{1S) = 12.1 ± 0.0 ± 0.2 ± 0.4 zs, (9.13) 

t{2S) = 21.6 ± 0.1 ± 0.4 ± 0.9 zs, and (9.14) 

r{3S) = 35.4 ± 0.4 ± 0.6 ± 1.7 zs, (9.15) 

where 1 zs (zeptosecond) is 10~^^ s. To the degree that the T system is non- 
relativistic, and therefore Equation 1.1 holds, the value of the bb wavefunction at 
the origin is 

|V'(0, 0,0)1^(15) = 53.0 ± 0.2 ± 0.8 fm-3, (9.16) 

1^(0,0,0)^(25) = 27.2 ± 0.2 ± 0.4 fm-3, and (9.17) 

1^(0,0,0)1^(35) = 20.9 ± 0.2 ± 0.3 fm-3. (9.18) 

The inverse cube roots of these values (0.27 fm, 0.33 fm, and 0.36 fm) characterize 
the physical sizes of the T mesons, though the form of the potential is needed to 
express them as the RMS of the wavefunction. These bottomonium |-?/'(0, 0, 0)p 
values, with ratios of 100:51:39, fall less rapidly from IS* to 35* than positronium, 
in which |-?/'(0, 0, 0)p have ratios of 100:35:19. We see in this the effect of the 
linear part of the Strong force potential: the T(3S') has more probability density 
at the origin than it would if its potential were purely Coulombic. In the sense 
that Fee is sensitive to the physical size and shape of the T meson, we have truly 
used a particle accelerator as a "giant microscope" to observe a structured object 
a quadrillion times smaller than a meter. 
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